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Abstract 

We address the decay rates of the energy for the damped wave equation when the damping 
coefficient b does not satisfy the Geometric Control Condition (GCC). First, we give a link with 
the controllability of the associated Schrodinger equation. We prove in an abstract setting that 
the observability of the Schrodinger group implies that the semigroup associated to the damped 
wave equation decays at rate 1/Vt (which is a stronger rate than the general logarithmic one 
predicted by the Lebeau Theorem). 

Second, we focus on the 2-dimensional torus. We prove that the best decay one can expect 
is 1/t, as soon as the damping region does not satisfy GCC. Conversely, for smooth damping 
coefficients b, we show that the semigroup decays at rate l/t^~^, for all e > 0. The proof relies 
on a second microlocalization around trapped directions, and resolvent estimates. 

In the case where the damping coefficient is a characteristic function of a strip (hence dis- 
continuous), Stephane Nonnenmacher computes in an appendix part of the spectrum of the 
associated damped wave operator, proving that the semigroup cannot decay faster than l/t'^^^. 
In particular, our study shows that the decay rate highly depends on the way b vanishes. 
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Part I 

The damped wave equation 

1 Decay of energy: a survey of existing results 

Let (M, g) be a smooth compact connected Riemannian d-dimensional manifold with or without 
boundary dM. We denote by A the (non-positive) Laplace-Beltrami operator on M for the metric 
g. Given a bounded nonnegative function, b € L°°{M), b{x) > on M, we want to understand the 
asymptotic behaviour as t — > +oo of the solution u of the problem 

dfu - Au + b{x)dtu = in R+ x M, 

u = on R+ X dM (if dM ^ 0), (1.1) 

(w, dtu)\t=i) = (uo, "i) in M. 

The energy of a solution is defined by 

E{u,t) = i(||V^i(t)||i.(,,) + ||a,u(t)||i.(,,,)). (1.2) 
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Multiplying (1.1) by dtu and integrating on M yields the following dissipation identity 

^E{u,t)^- f b\dtu\^dx, 
at Jj^i 

which, as b is nonnegative, implies a decay of the energy. As soon as 6 > C > on a nonempty open 
subset of M, the decay is strict and E{u, — >■ as i — >■ +00. The question is then to know at which 
rate the energy goes to zero. 

The first interesting issue concerns uniform stabilization: under which condition does there exist 
a function F{t), F{t) — > 0, such that 

E{u,t) < F{t)E{u,0) ? (1.3) 

The answer was given by Ranch and Taylor [RT74] in the case dM = and by Bardos, Lebeau and 
Ranch [BLR92] in the general case (see also [BG97] for the necessity of this condition): assuming 
that b S "^"(M), uniform stabilisation occurs if and only if the set {b > 0} satisfies the Geometric 
Control Condition (GCC). Recall that a set C M is said to satisfy GCC if there exists io > 
such that every geodesic 7 (resp. generalised geodesic in the case dM 7^ 0) of M with length larger 
than Lq satisfies 7 n w 7^ 0. Under this condition, one can take F{t) = Ce""* (for some constants 
C,K > 0) in (1.3), and the energy decays exponentially. Finally, Lebeau gives in [Leb96] the explicit 
(and optimal) value of the best decay rate k in terms of the spectral abscissa of the generator of the 
semigroup and the mean value of the function b along the rays of geometrical optics. 

In the case where {b > 0} docs not satisfy GCC, i.e. in the presence of "trapped rays" that do 
not meet {6 > 0}, what can be said about the decay rate of the energy? As soon as 6 > C > on 
a nonempty open subset of M, Lebeau shows in [Leb96] that the energy (of smoother initial data) 
goes at least logarithmically to zero (see also [Bur98]): 

E{u,t) < C{f{t)f (\\ua\\]j.(M)nHim + \Wi\?Hi(M)) ' for all t > 0, (1.4) 

with f{t) = j^^^ (where H^{M)nH^{M) and H^{M) have to be replaced by H^{M) and H^{M) 

respectively if dM = 0). Note that here, (/(t))^ characterizes the decay of the energy, whereas fit) 
is that of the associated semigroup. Moreover, the author constructed a series of explicit examples 
of geometries for which this rate is optimal, including for instance the case where M = is the 
two-dimensional sphere and {b > 0} D = 0, where N^^ is a neighbourhood of an equator of S^. 
This result is generalised in [LR97] for a wave equation damped on a (small) part of the boundary. 
In this paper, the authors also make the following comment about the result they obtain: 

"Notons toutefois qu'une etude plus approfondie de la localisation spectrale et des taux de 
decroissance de I'energie pour des donnees regulicres doit faire intervenir la dynamique globale du 
flot geodesique generalise sur M. Les theoremes [LR97, Theoreme 1] et [LR97, Theoreme 2] ne four- 
nissent done que les bornes a priori qu'on peut obtenir sans aucune hypothese sur la dynamique, en 
n'utilisant que les inegalites de Carleman qui traduisent "I'effet tunnel" ." 

In all examples where the optimal decay rate is logarithmic, the trapped ray is a stable trajectory 
from the point of view of the dynamics of the geodesic flow. This means basically that an important 
amount of the energy can stay concentrated, for a long time, in a neighbourhood of the trapped ray, 
i.e. away from the damping region. 

If the trapped trajectories are less stable, or unstable, one can expect to obtain an intermediate 
decay rate, between exponential and logarithmic. We shall say that the energy decays at rate f{t) 
if (1.4) is satisfied (more generally, see Definition 2.2 below in the abstract setting). This problem 
has already been adressed and, in some particular geometries, several different behaviours have been 
exhibited. Two main directions have been investigated. 

On the one hand, Liu and Rao considered in [LR05] the case where M is a square and the set 
{b > 0} contains a vertical strip. In this situation, the trapped trajectories consist in a family 
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of parallel vertical geodesies; these are unstable, in the sense that nearby geodesies diverge at a 
linear rate. They proved that the energy decays at rate ^ (i.e., that (1.4) is satisfied with 

f{t) = (^^) This was extended by Burq and Hitrik [BH07] (see also [Nis09]) to the case of 

partially rectangular two-dimensional domains, if the set {b > 0} contains a neighbourhood of the 
non-rectangular part. In [Phu07], Phung proved a decay at rate for some (unprecised) 5 > in 
a three-dimensional domain having two parallel faces. In all these situations, the only obstruction 
to GCC is due to a "cylinder of periodic orbits" . The geometry is flat and the unstabilities of the 
geodesic flow around the trapped rays are relatively weak (geodesies diverge at a linear rate). 

In [BH07] , the authors argue that the optimal decay in their geometry should be of the form j^W , 
for all e > 0. They provide conditions on the damping coefficient b{x) under which one can obtain 
such decay rates, and wonder whether this is true in general. Our main theorem (see Theorem 2.6 
below) extends these results to more general damping functions b on the two-dimensional torus. 

On the other hand, Christianson [ChrlO] proved that the energy decays at rate e~^^ for some 
C > 0, in the case where the trapped set is a hyperbolic closed geodesic. Schenck [Schll] proved an 
energy decay at rate e"*^* on manifolds with negative sectional curvature, if the trapped set is "small 
enough" in terms of topological pressure (for instance, a small neighbourhood of a closed geodesic), 
and if the damping is "large enough" (that is, starting from a damping function b, f3b will work for 
any /3 > sufficiently large). In these two papers, the geodesic flow near the trapped set enjoys 
strong instability properties: the flow on the trapped set is uniformly hyperbolic, in particular all 
trajectories are exponentially unstable. 

These eases confirm the idea that the decay rates of the energy strongly depends on the stability 
of trapped trajectories. 

One may now want to compare these geometric situations to situations where the Schrodinger 
group is observable (or, equivalently, controllable), i.e. for which there exist C > and T > such 
that, for all uq G L'^{M), we have 

holli.(M) <C TwVb e-^'^uo\\h^,,^dt. (1.5) 
Jo 

The conditions under which this property holds are also known to be related to stability of the 
geodesic flow. In particular, the works [BLR92], [LR05], [BH07, Nis09] and [ChrlO, Schll] can be 
seen as counterparts for damped wave equations of the articles [Leb92], [IIar89a, Jaf90], [BZ04] and 
[ARIO], respectively, in the context of observation of the Schrodinger group. 

Our main results are twofold. First, wc clarify (in an abstract setting) the link between the ob- 
servability (or the controllability) of the Schrodinger equation and polynomial decay for the damped 
wave equation. This follows the spirit of [Har89b] , [Mil05] , exploring the links between the different 
equations and their control properties (e.g. observability, controllability, stabilization...). More pre- 
cisely, we prove that the controllability of the Schrodinger equation implies a polynomial decay at 
rate for the damped wave equation (Theorem 2.3). 

Second, we study precisely the damped wave equation on the flat torus in case GCC fails. 
We give the following a priori lower bound on the decay rate, revisiting the argument of [BH07]: 
(1.1) is not stable at a better rate than j, provided that GCC is not satisflcd. In this situation, the 
Schrodinger group is known to be controllable (see [Jaf90] , [Kom92] and the more recent works [AMll] 
and [BZll]). Thus, one cannot hope to have a decay better than polynomial in our previous result, 
i.e. under the mere assumption that the Schrodinger flow is observable. 

The remainder of the paper is devoted to studying the gap between the a priori lower and upper 
bounds given respectively by j and on flat tori. For smooth nonvanishing damping coefficient 

b{x), we prove that the energy decays at rate ^4-? for all e > 0. This result holds without making 
any dynamical assumption on the damping coefficient, but only on the order of vanishing of b. It 
generalises a result of [BH07], which holds in the case where b is invariant in one direction. Our 
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analysis is, again, inspired by the recent microlocal approach proposed in [AMll] and [BZll] for 
the observability of the Schrodinger group. More precisely, we follow here several ideas and tools 
introduced in [MaclO] and [AMll]. 

In the situation where 6 is a characteristic function of a vertical strip of the torus (hence discon- 
tinuous), Stephane Nonnenmacher proves in Appendix B that the decay rate cannot be faster than 

. This is done by explicitly computing the high frequency eigenvalues of the damped wave oper- 
ator which arc closest to the imaginary axis (see for instance the figures in [AL03, AL12]). The fact 
that the decay rate 1 /t is not achieved in this situation was observed in the numerical computations 
presented in [AL12]. 

In contrast to the control problem for the Schodinger equation, this result shows that the stabi- 
lization of the wave equation is not only sensitive to the global properties of the geodesic flow, but 
also to the rate at which the damping function vanishes. 

2 Main results of the paper 

Our first result can be stated in a general abstract setting that we now introduce. We come back to 
the case of the torus afterwards. 

2.1 The damped wave equation in an abstract setting 

Let H and Y be two Hilbert spaces (resp. the state space and the observation/control space) with 
norms || • \\h and || • and associated inner products (•,•)-?/ ^-nd (•,-)i'- 

We denote by A : D{A) d H ^ H a nonnegative selfadjoint operator with compact resolvent, and 
B £ C{Y;H) a control operator. We recall that B* e £{H;Y) is defined by {B*h,y)Y = {h,By)H 
for all h E H and y G Y. 

Definition 2.1. We say that the system 

dtu + iAu^O, y^B*u, (2.1) 
is observable in time T if there exists a constant Kt > such that, for all solution of (2.1), we have 

ll"(0)||l^ < At r MtWydt. 
Jo 

We recall that the observability of (2.1) in time T is equivalent to the exact controllability in 
time T of the adjoint problem 

dtu + iAu^ Bf, w(0)=wo, (2.2) 

(see for instance [Leb92] or [RTTT05]). More precisely, given T > 0, the exact controllability in time 
T is the ability of finding for any uo,ui £ H a control function / G L^{0, T; Y) so that the solution 
of (2.2) satisfies u{T) = m. 

We equip H = D{Ai) x H with the graph norm 

\\{uo,u,)rn = \\{A + ld)iuo\\j, + \\u,\\j,, 

and define the seminorm 

l("o,wi)|« = WA^oWl + WuifH. 

Of course, if A is coercive on H, \ ■ j-^ is a norm on H equivalent to || • W-u- 

We also introduce in this abstract setting the damped wave equation on the space Ji, 

dfu + Au + BB*dtu^O, ^^^^ 
(m, dtu)\t=a = (uo, ui) e H, 
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which can be recast on first order system 



{ Z:^L.u.). "-{si)' ^-{"a -bb- ) . ^(-^) - ^(^) ^ ^<^*)^ (^^^) 

The compact injections D{A) ^ D{A^) ^ H imply that D{A) ^ Ti compactly, and that the 
operator A has a compact resolvent. 

We define the energy of solutions of (2.3) by 

E{u,t) = liWA^Wj, + WdtuWj,) = l\{u,dtu)\l.. 

Definition 2.2. Let / be a function such that /(t) when t +oo. We say that System (2.3) 
is stable at rate f{t) if there exists a constant C > such that for all (uq, ui) £ D{A), we have 

E{u,t)^ < Cf{t)\A{uo,ui)\'H, for all t > 0. 

If it is the case, for all k > 0, there exists a constant Cfc > such that for all {uq, ui) e D{A''), we 
have (see for instance [BD08, page 767]) 

E{u,t)i < Ck{f{t))''\\A''{uo,ui)\\n, foralH>0. 

Theorem 2.3. Suppose that there exists T > such that System (2.1) is observable in time T . Then 
System (2.3) is stable at rate 

Note that the gain of the log(t)5 with respect to [LR05, BH07] is not essential in our work. It is 
due to the optimal characterization of polynomially decaying semigroups obtained by Borichcv and 
Tomilov [BTIO]. 

This Theorem may be compared with the works (both presented in a similar abstract setting) 
[IIar89b] by Haraux, proving that the controllability of wave-type equations in some time is equivalent 
to uniform stabilization of (2.3), and [Mil05] by Miller, showing that the controllability of wave-type 
equations in some time implies the controllability of Schrodinger-type equations in any time. 

Note that the link between this abstract setting and that of Problem (1.1) is H = Y = L'^{M), 
A = -A with D{A) = H^{M) if DM = and H^{M) n i?o (M) otherwise, B is the muhiplication in 
L'^{M) by the bounded function y/b. 

As a first application of Theorem 2.3 we obtain a different proof of the polynomial decay results 
for wave equations of [LR05] and [BII07] as consequences of the associated control results for the 
Schrodinger equation of [Har89a] and [BZ04] respectively. 

Moreover, Theorem 2.3 provides also several new stability results for System (1.1) in particu- 
lar geometric situations; namely, in all following situations, the Schrodinger group is proved to be 
observable, and Theorem 2.3 gives the polynomial stability at rate ^ for (1.1): 

• For any nonvanishing b{x) > in the 2-dimensional square (resp. torus), as a consequence of 
[Jaf90] (resp. [MaclO, BZll]); for any nonvanishing b{x) > in the d-dimensional rectangle 
(resp. d-dimensional torus) as a consequence of [Kom92] (resp. [AMll]); 

• If M is the Bunimovich stadium and b{x) > on the neighbourhood of one half disc and on 
one point of the opposite side, as a consequence of [BZ04]; 

• If AI is a d-dimensional manifold of constant negative curvature and the set of trapped tra- 
jectories (as a subset of S*M, see [ARIO, Theorem 2.5] for a precise definition) has Hausdorff 
dimension lower than d, as a consequence of [ARIO]; 
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Moreover, Lebeau gives in [Leb96, Theoreme 1 (ii)] several 2-dimensional examples for which the 
decay rate jpg(2+t) optimal. For all these geometrical situations, Theorem 2.3 implies that the 
Schrodinger group is not observable. 

The proof of Theorem 2.3 relies on the following characterization of polynomial decay for Sys- 
tem (2.3). For z e C, we define on H the operator P{z) = A + Id+zBB* , with domain 
D{P{z)) ~ D{A). We prove in Lemma 4.2 below that P{is) is invertiblc for all s S M, s 7^ 0. 

Proposition 2.4. Suppose that 

for any eigenvector Lp of A, we have B*ip ^ 0. (2-5) 

Then, for all a > 0, the five following assertions are equivalent: 

The system (2.3) is stable at rate — , (2-6) 

There exist C >0 and so>0 such that for all s e R, \s\ > sq, \\{isld-A)^'^\\cCH) < C\s\i , (2.7) 
There exist C > and Sq > such that for all z Cz C, satisfying \z\ > Sq, 

and\Rc{z)\ < we have \\{zld-A)-^\\c(n) < C|Im(z)|^, (^.8) 

C| Im(z)| a 

There exist C > and sq > such that for all s £ R, |s| > sq, \\P{isy^\\c{H) < C|s|°"\ (2.9) 
There exists C > and Sq > such that for all s G M, |s| > Sq u € T^iA), 

wur^ < c{\s\i-^Pi^s)urH + \s\i\\B*ury). (2-^^) 

This proposition is proved as a consequence of the characterization of polynomial decay for general 
semigroups in terms of resolvent estimates given in [BTIO], providing the equivalence between (2.6) 
and (2.7). See also [BD08] for general decay rates in Banach spaces. Note in particular that the proof 
of a decay rate is reduced to the proof of a resolvent estimate on the imaginary axes. By the way, 
this estimate implies the existence of a "spectral gap" between the spectrum of A and the imaginary 
axis, given by (2.8). 

Note finally that the estimates (2.7), (2.9) and (2.10) can be equivalcntly restricted to s > 0, 
since P{—is)u = P{is)u. 



2.2 Decay rates for the damped wave equation on the torus 

The main results of this article deal with the decay rate for Problem (1.1) on the torus := 
(M/27rZ)^. In this setting, as well as in the abstract setting, we shall write P{z) — —A + z'^ + zb{x). 

First, we give an a priori lower bound for the decay rate of the damped wave equation, on 
T^, when GCC is "strongly violated", i.e. assuming that supp(6) does not satisfy GCC (instead of 
{b > 0}). This theorem is proved by constructing explicit quasimodes for the operator P{is). 

Theorem 2.5. Suppose that there exists (xo,^o) G T*T^, ^0 7^ 0, such that 

{b > 0} n {xo + rCo , T e K} = 0. 

Then there exist two constants C > and kq > such that for all n G N, 

ll^(*^^0)"'||£(L^(T^)) > C. (2.11) 

As a consequence of Proposition 2.4, polynomial stabilization at rate for e > is not 
possible if there is a strongly trapped ray (i.e. that does not intersect supp(6)). More precisely, 
in such geometry. Theorem 2.5 combined with Lemma 4.6 and [BD08, Proposition 1.3] shows that 
"^i(^) ^ i+t J for some C > (with the notation of [BD08] where mi{t) denotes the best decay rate). 

Then, the main goal of this paper is to explore the gap between the a priori upper bound for 

the decay rate, given by Theorem 2.3, and the a priori lower bound j of Theorem 2.5. Our results 
are twofold (somehow in two opposite directions) and concern either the case of smooth damping 
functions b, or the case b = lu, with t/ C T^. 



7 



2.2.1 The case of smooth damping coefficients 

Our main result deals with the case of smooth damping coefficients. Without any geometric assump- 
tion, but with an additional hypothesis on the order of vanishing of the damping function b, we prove 
a weak converse of Theorem 2.5. 

Theorem 2.6. Let M = with the standard flat metric. There exists Eq > satisfying the following 
property. Suppose that b is a nonnegative nonvanishing function on satisfying \fb e "^""(T^) and 
that there exist e g (0, Eq) and > such that 

\Vb{x)\ < Ceb^-%x), for x e T^. (2.12) 

Then, there exist C > and sq > such that for all s G R, \s\ > sq, 

\\P{is)~'^\\ciLHT^)) <C\s\\ withS = 8E (2.13) 

As a consequence of Proposition 2.4, in this situation, the damped wave equation (1.1) is stable at 
rate — \— ■ 

Following carefully the steps of the proof, one sees that Eq = ^ works, but the proof is not 
optimized with respect to this parameter, and it is likely that it could be much improved. 

One of the main difficulties in understanding the decay rates is that there exists no general 
monotonicity property of the type "bi{x) < b2{x) for all x =^ the decay rate associated to the 
damping 62 is larger (or smaller) than the decay rate associated to the damping 61". This makes a 
significant difference with observability or controllability problems of the type (1.5). 

Assumption (2.12) is only a local assumption in a neighbourhood of d{b > 0} (even if it is stated 
here globally on T^). Far from this set, i.e. on each compact set {b > bo} for bo > 0, the constant 
Ce can be choosen uniformly, depending only on bo, and not on e. Hence, e somehow quantifies the 
vanishing rate of the damping function b. 

An interesting situation is when the smooth function b vanishes like e~^ in smooth local coordi- 
nates, for some a > 0. In this case. Assumption (2.12) is satisfied for any e > 0, and the associated 
damped wave equation (1.1) is stable at rate for any 6 > 0. This shows that the lower bound 
given by Theorem 2.5, as well as the decay rate j, are sharp in general. This phenomenon had 
already been remarked by Burq and Hitrik in [BH07] in the case where b is invariant in one direction. 

Typical smooth functions not satisfying Assumption (2.12) are for instance functions vanishing 
like sin(i)^e~i. We do not have any idea of the decay rate achieved in this case (except for the a 
priori bounds and i). 

Theorem 2.6 generalises the result of [BH07], which only holds if b is assumed to be invariant 
in one direction. Our proof is based on ideas and tools developped in [MaclO, AMll] and especially 
on two-microlocal semiclassical measures. One of the key technical points appears in Section 13: 
we have to construct, for each trapped direction, a cutoff function invariant in that direction and 
adapted to the damping coefficient b. We do not know how to adapt this technical construction 
to tori of higher dimension, d > 2; hence we do not know whether Theorem 2.6 holds in higher 
dimension (although we have no reason to suspect it should not hold). Only in the particular case 
where b is invariant in d — 1 directions can our methods (or those of [BH07]) be applied to prove the 
analogue of Theorem 2.6. 

Note that if GCC is satisfied, one has (on a general compact manifold M) for some C > 1 and 
all |s| > So the estimate 

\\Pits)-'\\ciLHM))<C\s\-\ (2.14) 

instead of (2.13). Estimate (2.14) is in turn equivalent to uniform stabilization (see [Hua85] together 
with Lemma 4.6 below). 
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Remark 2.7. As a consequence of Theorem 2.6 on the torus, we can deduce that the decay rate 
t^T+^ also holds for Equation (1.1) if M = (0, tt)^ is the square, one takes with Dirichlet or Neumann 
boundary conditions, and the damping function b is smooth, vanishes near dM and satisfies As- 
sumption (2.12). First, we extend the function b as an even (with respect to both variables) smooth 
function on the larger square (— 7r,7r)^, and using the injection i : (— 7r,7r)^ — T^, as a smooth func- 
tion on T^, still satisfying (2.12). Moreover, D{Ad) (resp. D{An)) on (0,7r)^ can be identified as 
the closed subspace of odd (resp. even) functions of D{Ad) (resp. D(Ajv)) on (— tt, tt)^. Using again 
the injection i, it can also be identified with a closed subspace of i7^(T^). The estimate 

lk|U2(T2) < C|sh|P(is)M|U2(T2) for all u e H^{T^), 

is thus also true on the square (0,7r)^ for Dirichlet or Neumann boundary conditions. In particular, 
this strongly improves the results of [LR05]. 

The lower bound of Theorem 2.5 can be similarly extended to the case of a square with Dirichlet 
or Neumann boundary conditions, implying that the rate j is optimal if GCC is strongly violated. 

2.2.2 The case of discontinuous damping functions 

Appendix B (by Stephane Nonnenmacher) deals with the case where b is the characteristic function 
of a vertical strip, i.e. b ~ Btu, for some B > Q and ?7 = (a, 6) x T C T^. Due to the invariancc of b 
in one direction, the spectrum of the damped wave operator A splits into countably many "branches" 
of eigenvalues. This structure of the spectrum is illustrated in the numerics of [AL03, AL12]. 

The branch closest to the imaginary axis is explicitly computed, it contains a sequence of eigen- 
values (zi)igN such that Imz^ — > oo and |Rezi| < ^i^^°y/2 ■ This result is in agreement with the 
numerical tests given in [AL12]. 

As a consequence, for any e > and C > 0, the strip ||Rcz| < C| Im(z)|~'^/^"'"'^} contains in- 
finitely many poles of the resolvent (z Id —A)~^, so item (2.8) in Proposition 2.4 implies the following 
obstruction to the stability of this damped system : 

Corollary 2.8. For any e > 0, the damped wave equation (1.1) on with the damping function 
(B.l) cannot be stable at the rate jr/kr^ ■ 

The same result holds on the square with Dirichlet or Neumann boundary conditions. 

More precisely, in this situation. Lemma 4.6 and [BD08, Proposition 1.3] yield that mi{t) > 
(14.^)2/3 J for some C > (with the notation of [BD08] where mi (t) denotes the best decay rate) . 

This corollary shows in particular that the regularity conditions in Theorem 2.6 cannot be com- 
pletely disposed of if one wants a stability at the rate 1/t^^'^ for small e. 

2.3 Some related open questions 

The various results obtained in this article lead to several open questions. 

1. In the case where b is the characteristic function of a vertical strip, our analysis shows that the 
best decay rate lies somewhere between and but the "true" decay rate is not yet clear. 

42 43 

2. It would also be interesting to investigate the spectrum and the decay rates for damping 
functions b invariant in one direction, but having a less singular behaviour than a characteristic 
function. In particular, is it possible to give a precise link between the vanishing rate of b and 
the decay rate? 

3. In the general setting of Section 2.1 (as well as in the case of the damped wave equation on the 
torus), is the a priori upper bound for the decay rate optimal? 

45 

4. For smooth damping functions vanishing like e~^, Theorem 2.6 yields stability at rate 

for all (5 > 0. Is the decay rate 7 reached in this situation? Can one find a damping function b 

4 ^ 

such that the decay rate is exactly j7 



9 



5. The lower bound of of Theorem 2.5 is still valid in higher dimensional tori. Is there an analogue 
of Theorem 2.6 (i.e. for general "smooth" damping functions) for T"^, with c? > 3? 

Part II 

Resolvent estimates and stabilization in 
the abstract setting 

3 Proof of Theorem 2.3 assuming Proposition 2.4 

To prove Theorem 2.3, we express the observability condition as a resolvent estimate (also known as 
the Hautus test), as introduced by Burq and Zworski [BZ04], and further developed by Miller [Mil05] 
and Ramdani, Takahashi, Tcnenbaum and Tucsnak [RTTT05]. For a survey of this notion, wc refer 
to the book [TW09, Section 6.6]. 

In particular [Mil05, Theorem 5.1] (or [TW09, Theorem 6.6.1]) yields that System (2.1) is ob- 
servable in some time T > if and only if there exists a constant C > such that we have 

< C(||(yl- Ald)u|||^ + \\B*u\\l-), for all A e K and u e D{A). 

As a first consequence. Assumption (2.5) is satisfied and Proposition 2.4 applies in this context. 
Moreover, we have, for all s S K and u <S D{A), 

\\u\\h < C{\\{A- Id +isBB* - isBB*)u\\'jj + \\B*u\\l.) 

< C{\\P{is)ufH + s^\\BB*ufH + I|B*w|Iy) (3.1) 
Since B <E C{Y; H), we obtain for s > 1 and for some C > 0, 

< C{\\P{is)u\\l + s^B*ury) < C{s^P{is)u\\l + s^B*u\\l.). 

Proposition 2.4 then yields the polynomial stability at rate ^ for (2.3). This concludes the proof of 
Theorem 2.3. □ 

4 Proof of Proposition 2.4 

Our proof strongly relies on the characterization of polynomially stable semigroups, given in [BTIO, 
Theorem 2.4], which can be reformulated as follows. 

Theorem 4.1 ([BTIO], Theorem 2.4). Let (e*-^)t>o be a bounded -semigroup on a Hilbert space 
%, generated by A. Suppose that H Sp(^) = 0. Then, the following conditions are equivalent: 

\\e'^A-^\\c(H)=0[t-°'), ast^+oo, (4.1) 

\\{isld-A)-^\\^(^^^^Oi\s\-^), ass-^oo. (4.2) 

Let us first describe some spectral properties of the operator A defined in (2.4). 
Lemma 4.2. The spectrum of A contains only isolated eigenvalues and we have 

Sp(^)c + u ([-11^11 + 

with keiiA) = ker(A) x {0}. 
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Moreover, the operator P{z) is an isomorphism from D{A) onto H if and only if z ^ Sp(,4). // 
this is satisfied, we have 

^^^"^ ^ [ Piz)-^izBB* +z^ld)-Id zP{z)-^ )■ ^^■'^> 

The localization properties for the spectrum of A, stated in the first part of this lemma are 
illustrated for instance in [AL03] or [AL12]. 

This Lemma leads us to introduce the spectral projector of A on ker(^), given by 



^0 = 7^ [ {zid-Ay^dz e c{n), 



where 7 denotes a positively oriented circle centered on with a radius so small that is the single 
eigenvalue of A in the interior of 7. We set H = (Id — no)'H and equip this space with the norm 

||(uo,"i)||^ := IK,Ui)l?i = U^uoWh + WMh, 

and associated inner product. This is indeed a norm on H since ||(mo,wi)||^ = is equivalent to 
(wo,Mi) G ker(A) x {0} = noH. 

Besides, we set A ~ A\^ with domain D{A) = D{A) OH. A first remark is that Sp(^) = 
Sp(yl) \ {0}, so that Sp(y4) n = 0. 

The remainder of the proof consists in applying Theorem 4.1 to the operator A in H. We first 
check the assumptions of Theorem 4.1 and describe the solutions of the evolution problem (2.4) (or 
equivalently (2.3)). 

Lemma 4.3. The operator A generates a contraction '^^'^ -semigroup on ik., denoted (e*-^)t>o. More- 
over, for all initial data Uq £ Ji, Problem (2.4) (or equivalently (2.3)j has a unique solution 
U € {W^ ; T-L) , issued from Uq, that can he decomposed as 

U(t) = e*-^(Id-no)?7o + HoC/o, for all t > 0. (4.4) 

As a consequence, we can apply Theorem 4.1 to the semigroup generated by A. The proof of 
Proposition 2.4 will be achieved when the following lemmata are proved. 

Lemma 4.4. Conditions (2.6) and (4.1) are equivalent. 

Lemma 4.5. Conditions (2.9) and (2.10) are equivalent. Conditions (2.7) and (2.8) are equivalent. 
Lemma 4.6. There exist C > 1 and sq > such that for s e |s| > sq, 

|i(zsId-^)-i||^(^) - 1^ < \\izsld-A)-'\\cm < ||(*sld-i)-i||^(^) + (4.5) 

and 

C-'\s\\\P{is)-^ciH) < Uisld^A)-^cm) <C{1 + \s\\\P{is)-^\\ciH)) ■ (4.6) 

In particular this implies that (4.2), (2.7) and (2.9) are equivalent. 

The proof of Lemma 4.6 is more or less classical and we follow [Lcb96, BH07]. 

Proof of Lemma 4-2. As A has compact resolvent, its spectrum contains only isolated eigenvalues. 
Suppose that z E Sp(^); then we have, for some (mo,wi) £ D{A) \ {0}, 

Ul = ZUq, 

-Auq — BB*ui = zui, 
and in particular 

Auo + z^uq + zBB*uo = 0, (4.7) 
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with Mo e D{A) \ {0}. 

Suppose that z G iM., then, this yields Auq — lm{z)^uo + ilm{z)BB*uo = 0. Following [Leb96], 
taking the inner product of this equation with uq yields iIm(z)||i?*Mo||y = 0. Hence, either Im(z) = 0, 
or B*uq = 0. In the first case, Auq = 0, i.e. uq G ker(A), and ui = 0. This yields ker(^) C 
ker(A) X {0} (and the other inclusion is clear). In the second case, uo is an eigenvector of A associated 
to the eigenvalue lm{z)^ and satisfies B*uo = 0, which is absurd, according to Assumption (2.5). 
Thus, Sp(^) n « C {0}. 

Now, for a general eigenvalue z G C, taking the inner product of (4.7) with uq yields 



(4.8) 



{Auo, uo)h + (Re(z)2 - Im(z)2)||uo|||j + Reiz)\\B*uo\\^Y = 0, 
2Rc{z)lm{z)\\uQ\\jf + Im(z)||B*Mo|lY = 0. 

If Im(z) 0, then, the second equation of (4.8) together with Sp(yl) n iM C {0} gives 

If Im(z) = 0, then, the first equation of (4.8) together with {Auq, uo)h > gives — Re(z)||i?*ito||y > 
Re(2;)2||Mo||l/, which yields 

0>Re{z)>-\\B*\\l^H;Yy 

Following [Leb96], we now give the link between P{z)~^ and (zld— ^)~^ for z ^ Sp(^). Taking 
F ~ (/o, /i) G "H, and U = {uo,ui), we have 

(z Id -yl) (7^ ( pi = (4.9) 
[ P{z)uo = fi + [BB* + zld)fo. 

As a consequence, we obtain that P{z) : D{A) — > is invertible if and only if (zld— .4) : D{A) — > H 
is invertible, i.e. if and only if z ^ Sp(.4). Moreover, for such values of z, System (4.9) is equivalent 
to 

Mo = P{z)-^fi + P{z)-\BB* + zId)/o, 
Ml = zP{z)-\fi + zP{z)-\BB* + zId)/o - /o, 

which can be rewritten as (4.3). This concludes the proof of Lemma 4.2. □ 

Proof of Lemma 4-3. Let us check that ^ is a maximal dissipative operator on H [Paz83]. First, it 
is dissipative since, for U = (mo,mi) e D{A), 

{AU,U)^ = iA^ui,A^uo)H - {Auo,ui)h - {BB*ui,ui)h = -\\B*ui\\l < 0. 

Next, the fact that ^ — Id is onto is a consequence of Lemma 4.2. Hence, for all F E T-l C Ti, 
there exists U <E D{A) such that {A — ld)U = F. Applying (Id— Ho) to this identity yields (A — 
Id)(Id — Ho)?/ = F, so that yl— Id : D{A) — "H is onto. According to the Lumer-Phillips Theorem (see 
for instance [Paz83, Chapter 1, Theorem 4.3]) A generates a contraction "^"-semigroup on H. Then, 
Formula (4.4) directly comes from the linearity of Equation (2.4) (or equivalently (2.3)) together 
with the decomposition of the initial condition J/o = (-^ ^ no)[/o + Ho?7o- D 

Proof of Lemma 4-4- Condition (4.1) is equivalent to the existence of C > such that for all t > 0, 
and Uo & Ti, we have 



We'-^A-'UoWn < ^IIC/oIIh. 

This can be rephrased as 



We'^ifoWn < ^Mf/olk, (4.10) 

for all t > 0, and Uq G D{A). Now, take any Uq — (mo,mi) e D{A), and associated projection 
ilo = (Id-no)?7o e D{A). According to (4.4), we have 

E{u,t) = liWA^uml + Wdtuiml) = he'^Uo + HoC/o|?, = ke'^UoWl,, 
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and 

\AUo\n = lAUo + AUoUoln = WAUoWn- 
This shows that (4.10) is equivalent to (2.6), and concludes the proof of Lemma 4.4. 



□ 



Proof of Lemma 4-5. First, (2.9) clearly implies (2.10). To prove the converse, for u e D{A), we 
have 



(P(is)it, u)h = {{A - Id)u, u)jj+ is\\B*u\\ 



Taking the imaginary part of this identity gives s||i3*uj|y = Im(P(is)w, u)h; so that, using the Young 
inequality, we obtain for all £ > 0, 



\s\nB*u\\l = \s\-^-' 



\Iiti{P{is)u,u)h\ < 



4£ 



-\ms)urH+e\\u''' 



H- 



Plugging this into (2.10) and taking e sufficiently small, we obtain that for some C > and so > 0, 
for any s G R, \s\ > sq, 

\\u\\l<C\s\i-^P{is)u\\l, 
which yields (2.9). Hence (2.9) and (2.10) are equivalent. 

Second, Condition (2.8) clearly implies (2.7) and it only remains to prove the converse. For 2; e C, 
we write r = Re(z) and s = Im(z). We have the identity 



Hence, assuming 
this gives 



{{r + is)ld-A)-^ = {isld-A)-\ld+r{isld-A)-^) \ 



\\r{isld-Ar'\\c(n) < 



(4.11) 
(4.12) 



'ld+r{isld-Ay 



C{H) 



Y,{-l)\r{^sld-Ar'] 



k=0 



< 2. 



As a consequence of (4.11) and (2.7), we then obtain 

\\iir + is)ld-A)-'\\^^^^ < 2\\iisld-A)-'\\cin) < 2C|s|i, 
for all s > sq, under Condition (4.12). Finally, (2.7) also yields 

Wriisld-Ar'Wan) < 

so that Condition (4.12) is realised as soon as |r| < — This proves (2.8) and concludes the 
proof of Lemma 4.5. □ 



Proof of Lemma 4-6- To prove (4.5), we first remark that the norms 



and 



Ifi are equiv- 



alent on so that the norms 



£(«) 



and 11 • WccH) ^I'S equivalent on £('H). Next, we have 



iisld-A)-\ld-no) = {isld-A)-^{ld-na) and 

\\{2sld-A)-'\\c(n) = \\{tsld-Ar\ld-Ilo)\\c(n) = \\{tsId-A)-\ld -no)\\cin) 
< \\ii.sld~A)-'\\c(H) + Uisld~A)-'Uo\\ciH), 



together with 



Wizsld-Ar'Wan) = !l(*.sId-i)"i(Id-Ho) + (*,sId-^)-iHo|U(„) 
< \\{isld-A)-'\\c(n) + \\iisld-A)-'no\\c(n)- 
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Moreover, for |s| > 1, we have 

\\{isld-A)-'no\\cin) = \\{is)-'llo\\ciH) = ^l|no||£(«) - j^, 
which concludes the proof of (4.5). 

Let us now prove (4.6). For concision, we set Hi = D{A'^) endowed with the graph norm 
= \\{A + Id)2u||_f/ and denote by iJ-i = D{A2y its dual space. The operator A can be 
uniquely extended as an operator £{Hi; H-i), still denoted A fo simplicity. With this notation, the 
space H_i can be equipped with the natural norm = ||(^ + ld)~^u\\ff. 

As a consequence of Formula (4.3), and using the fact that Sp(^)niR C {0}, there exist constants 
C > 1 and So > such that for all s S R, \s\ > sq, 

C-^M{s) < \\{tsld-A)-^cin) < CM{s) (4.13) 

with 

M{s) ={\\P{is)-HBB* +isld)\\ciH,) + \\P{^s)-^\\c(H-M,) 

+ \\P{is)-\isBB* - sHA) - Id \\ciHr,H) + \\sP{is)-^\\c(H)) (4.14) 
On the one hand, this dirccly yields for s E M, \s\ > sq, 

\s\\\Piisr'\\ciH) < C\\iisU-A)-'\\ciH)- 

This proves that (4.2) implies (2.9). 

On the other hand, we have to estimate each term of (4.14). First, using Au = P{is)u + s^u — 
isBB*u, we have 

II^^IIhi = ll^^^illlf + WA'h = {P{is)u + s^u - isBB*u,u)^ + \\u\\jj 

= Re {P{ts)u,u)^ + {s^ + l)\\u\\l < C {\\P{is)ufH + {s^ + l)\H\l) 
<c{l + is' + l)\\P{is)-YciH)) \\Pi^s)u\\j„ 

so that 

\\Pi^sr'\\c(H;H,) < C (1 + (|s| + l)\\P{t.s)-'\\ciH)) ■ (4.15) 

Second, the same computation for {P{is)^^)* = (A— Id—isBB*)^^ (the adjoint of P{is)^^ in 
the space H) in place of P(is)~^ leads to {P{is)^^)* e C{H; Hi), together with the estimate 

\m^sr'r\\ciH:H,) <c{i + {\s\ + mpiisywaH)) . 

By transposition, we have *(P(is)^^)* e C{H^i; H), together with the estimate 

\\HP{^s)-'r\\ciH^,,H) < \m^s)-'r\\c(H■M^) <C{1 + (|s| + l)\\P{zs)-'\\ciH)) . (4.16) 
Moreover, *(P(is)~^)* is defined, for every u G H,v G by 

{'{P{ts)-'rv,u)^ = {v, {P{^s)-'ru)^_^^^^ = [{A + Uyh, {A + ld)HP{^s)-'ru)^. 
In particular, taking v G H gives 

{\p{ts)-'rv,u)^ = {p{tsr\,iL)^, 

which implies that the restriction of the operator *(P(is)^^)* to H coincides with P{is)~^. For 
simplicity, we will denote P{is)~^ for *{P{is)~-^)* . 



14 



Equation (4.16) can thus be rewritten 

\\Pi^s)-'\\ciH^,.,H) < C (1 + (|.| + ^)\\Pi^s)-^ciH)) . (4.17) 
Then, we have P{is)~^{isBB* - Id) - Id = P{is)-^A, so that 

\\Pits)-\lsBB* - sHd) -IdWciHuH) = \\Piis)-'A\\c(H,-H) < \\Piis)-'\\c(H^uH)\\A\\c(HuH^,) 

<{l + i\s\ + l)\\Piis)-^c(H)) (4.18) 

Third, for |s| > 1 we write 

P(is)-\BB* + is Id) = - (P(is)~^A - Id) , (4.19) 

s 

and it remains to estimate the term \\P{is)^^A\\c(^H-^-j in (4.14). For / G Hi, we set u = P{is)^-^Af. 
We have u E Hi, together with 

{A-s^ ld+isBB*)u = Af. 
Taking the real part of the inner product of this identity with u, we find 

\\A-^u\\l - s^\\u\\j, = Re{Af,u)H < \\Af\\H_M\H, < C\\f\\HAW\\m, 
as A e C{Hi,H_i). Hence 

\\u\\l,<C{l + s')\\u\\l + C\\f\\%^ 

Using (4.17), this gives 

llt^ill^, < C(l + s^)\\Pi^s)-'A\\l^^^.^^^\\frH^ + CWfWl^ 

< C{1 + s')\\P{^s)-A\l^^_^.^^^\\frH^ + CWfrn^ 

< C(l + s') (1 + i\s\ + l)\\Pi^s)-A\ciH))' WfWl,, 

and finaUy WPiisy^ A\\c[Hi) < C(l + \s\) (l + (|s| + l)\\P{is)-A\c{H)) ■ Coming back to (4.19), we 
have, for \s\ > 1, 

\\Pits)-\BB* + isld)\\c(H,) <C{1 + \s\\\P{is)-A\ciH)) . (4.20) 
Finahy, combining (4.15), (4.18) and (4.20), together with (4.13)-(4.14), we obtain for |s| > 1, 

\\itsld-A)-'\\c(n) <C{1 + \s\\\P{zs)-^c(H)) ■ 
This concludes the proof of Lemma 4.6. 

□ 

Part III 

Proof of Theorem 2.6: smooth damping 
coefficients on the torus 

To prove Theorem 2.6, we shall instead prove Estimate (2.9) with a = (which, according to 
Proposition 2.4, is equivalent to the statement of Theorem 2.6). Let us first recast (2.9) with a = 
in the semiclassical setting : taking /i = s^^, we arc left to prove that there exist C > 1 and ho > 
such that for all h < ho, for all u G i7^(T^), we have 

\\u\\L2il2) < Ch-'\\Pii/h)u\\L2^j2) (4.21) 
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We prove this inequality by contradiction, using the notion of semiclassical measures. The idea 
of developing such a strategy for proving energy estimates, together with the associate technology, 
originates from Lebeau [Leb96]. 

We assume that (4.21) is not satisfied, and will obtain a contradiction at the end of Section 11. 
Hence, for all n g N, there exists < /i„ < ^ and m„ £ H^{T^) such that 

Tl 

l|Un||L2(T2) > -rg-\\P{i/hn)Un\\L^(T2). 

Setting v„ = u„/||u„||i2(T2), and 

Pi'" = -hlA - 1 + ih^bix) = hlP{i/K), 

we then have, as n — >■ oo, 

K 0+, 

!|Wn||L2(T2) = 1, 
^T'-'lln''"««||L^(T^)->0. 

Our goal is now to associate to the sequence (w„, hn) a semiclassical measure on the cotangent bundle 
on T*T2 = T2 X (M^)* (where (M^)* is the dual space of M^). To obtain a contradiction, we shah 
prove both that //(r*!^) = 1, and that ^ = on T*?^. 

From now on, we drop the subscript n of the sequences above, and write h in place of /i„ and vt 
in place of Vn- We study sequences (/i, Vh) such that /i — >■ 0+ and 

||w,,||i2(T2) = 1 

\\P^Vh\\mT^-)^o{h^+'), as/i^O+. ^' ' 

In particular, this last equation also yields the key information 

{hvh,Vh)L^{Y^) = lm{Pl^Vh,Vh) L^(T^) = o(/i^+'^'), as /i ^ 0+. 

In the following, it will be convenient to identify (M^)* and through the usual inner product. 
In particular, the cotangent bundle T*t^ = x (R^)* will be identified with x R^. 



5 Semiclassical measures 

We denote by T*T'^ the compactification of T*T'^ obtained by adding a point at infinity to each fiber 
(i.e., the set x (R^ U {oo})). A neighbourhood of (x,oo) e T*V is a set [/ x ({oo} U R^ \ iC), 
where t/ is a neighbourhood of x in and K a compact set in R^. Endowed with this topology, the 
set T*T2 is compact. 

We denote by S^{T*T'^), 5" for short, the space of functions a{x,^) that satisfy the following 
properties: 

1. a e <^#'°°(T*T2). 

2. There exists a compact set iiT C R^ and a constant fco G C such that a(a;, ^) = fco for all 

^ e r2 \ 

Note that we have in particular '^^{T*T'^) C 5'°(T*T2). 

To a symbol a £ 5°(T*T^), we associate its semiclassical Weyl quantization Op,j(a) by For- 
mula (A.l), which, according to the Calderon-Vaillancourt Theorem (see Appendix A) defines a 
uniformly bounded operator on L^(T^). 

From the sequence {vh , h) (see for instance [GL93] ) , we can define (using again the Calderon- 
Vaillancourt Theorem) the associated Wigner distribution € [S^)' by 

("^'''«)(s«)'.so = {'^Vh{o)vh,Vh)L2^^2) , for ah a £ S"{T*t'). (5.1) 
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Decomposing Vh and a in Fourier series, 

It: Jf2 Zn J12 
the expression (5.1) can be more explicitly rewritten as 

Proposition 5.1. The family (V^) is bounded in (S^)' . Hence, there exists a subsequence of the 

sequence {h,Vh) and an element fi £ (S^)' , such that fi weakly in (S^)' , i.e. 

(Op^(a)w/i,w,i)^2(T2) (m, a) (SO)', so " ^ 5°(T*T^). (5.2) 

In addition, {fi, a) ^goy go is nonnegative if a is; in other words, fi may be identified with a nonnegative 
Radon measure on T*T^. 

Notation: in what follows we shall denote by Ai~^ iT*T'^) the set of nonnegative Radon measures 
on T*f2. 

Proof. The proof is an adaptation from the original proof of Gerard [Gcr91] (see also [GL93] in the 
scmiclassical setting). 

The fact that the Wigner distributions V'^ are uniformly bounded in (5")' follows from the 
Calderon-Vaillancourt theorem (see Appendix A), and from the boundedness of (w/j) in L^(T^). 
The sharp Garding inequality gives the existence of C > such that, for all a > and h > 0, 

(Op,,(a)uh,w/0i2(T2) > -C/i||uh||i2(T2), 
so that the distribution fi is nonnegative (and is hence a measure). 

□ 



6 Zero-th and first order informations on 

To simplify the notation, we set 

Pi; = P^ + ihb{x), with P',' = -h'A - 1 = Op,(|e|' - 1). 

The geodesic flow on the torus : T*T'^ — >• T*T'^ for t £ M is the flow generated by the 
Hamiltonian vector field associated to the symbol ^(|CP — l)j i-c by the vector field ^ • 9^ on T*T^. 
Explicitely, we have 

MxA) = ix + T^,0. reR, ix,OeT*T\ 

Note that (j)r preserves the ^-component, and, in particular every energy layer = C > 0} C r*T^. 
Now, we describe the first properties of the measure /i implied by (4.22). 

We recall that for v e &'{T*T^), {(f)r)*i^ £ &'{T*T^) is defined by ((^^)*zy,a) = {v^aocj)^) for 
all a € '^^°°(T*T^). In particular, {(f>r)*v is a measure if i/ is. We shall say that u is an invariant 
measure if it is invariant by the geodesic flow, i.e. {(Pt)*^ = v for all r G K. 

Proposition 6.1. Let fi be as in Proposition 5.1. We have 

1. supp(/.j) C = 1} (hence is compact in T*T'^), 

2. fi{T*T^) = 1, 

3. fi is invariant by the geodesic flow, i.e. {(j)r)*lJi- — 

4-. (/^: ^);vi^(T*T2) <;g'0(T*T2) = 0, where M^T*^"^) denotes the space of compactly supported mea- 
sures on T*T^. 
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In other words, fi is an invariant probability measure on T*T'^ vanishing on {b > 0}. 

These are standard arguments, that we reproduce here for the reader's comfort. In particular, 
we recover aU informations required to prove the Bardos-Lebeau-Rauch- Taylor uniform stabihzation 
theorem under GCC. But we do not use here the second order informations of (4.22); this will be 
the key point to prove Theorem 2.6. 



Proof. First, we take x £ (T*T'^) depending only on the ^ variable, such that x > 0, xiO = 

x( 
W 



for 1^1 < 2, and xiO = 1 for |CI > 3. Hence, jfMj G 'if°°{T*T^) and we have the exact composition 



formula 



Op,(x) = Op, ( -^L) F/' 



' 



since both operators are Fourier multipliers. Moreover, Op;j ^ ^ bounded operator on 

i^(T^). As a consequence, we have 

together with 



Op,. (^) - ■'■ («^' (#Tj ,,,, 

Since ||i2(-x2) = o(l) and ||u/i||l2(t2) = 1, both terms in this expression vanish in the limit 

h — >■ 0+. This implies that (/i, x)m{ T'T^ } '€°(T*r^) ~ ^i^ice this holds for all x above, we have 
supp(//) C = 1}, which proves Item 1. 

In particular, this implies that ^, [t*V \ T*T^^ = 0. Now, Item 2 is a direct consequence of 1 = 
||i'ft,|li2(T2) (a*, I)a^(t^) 'ro(T^) s-'^d Item 1. Item 4 is a direct consequence of (6t;/i, w/i)^2(t2) = 

0(1). 

Finally, for a G ^^°°(T*T2), wc rccaU that 

[Po", Op, (a)] = - Op,({|CP - 1, a}) = - Opje • d^a), 

is a consequence of the Weyl quantization (any other quantization would have left an error term of 
order OQi^)). Hence, (5.1) yields 

together with 

(y'\e-a,a)^,(y.^2)_^oo(r.T2) = ^([-Po'',Op,(a)] vh^vh)^,^^,^ 

" 'y^^hia)vh,P^Vh)^^^^^^ - ^ (Op,(a)F(^'w/^,w/^)^2(T2) 

" ('^^hi<^^^h,PbVh)^,_^^^^ - ^ (Op,(a)Pb''uh,Uh)^2(T2) 

- \ (Op,(a)uh,&w/0^2(T2) - ^ (Op,(a)6v,,u,)^2(T2) ■ (6.2) 

In this expression, we have [Opfj {a)vh, Pj^Vh) j^2(j2-f ~^ and (Op,(a)P^''u/i, w/i)^2(t2) ~^ since 
\\PhVh\\L^{i:^) = o{h). Moreover, the last two terms can be estimated by 

I {OPhiO')bvh,Vh)L2^j2~i I < II V6w/i||L2(T2)||V50p,j(a)w/i||L2(T2) = o(l), (6.3) 
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since '!;ft)^2(T2) = o(l). This yields {V^ ■ d^a) 0, so that, using (6.1), 
(^,^ • dxo) j^(j^,f2-^ <g'0(T'T2) = for all a G '^^°°(T*T^). Replacing a hy a o (p^. and integrating with 
respect to the parameter r gives {4)t)*H = which concludes the proof of Item 3. 

□ 



7 Geometry on the torus and decomposition of invariant mea- 
sures 

7.1 Resonant and non-resonant vectors on the torus 

In this section, we collect several facts concerning the geometry of T*T'^ and its resonant subspaces. 
Most of the setting and the notation comes from [AMll, Section 2]. 

We shall say that a submodule A C is primitive if (A) n = A, where (A) denotes the linear 
subspace of spanned by A. The family of all primitive submodules will be denoted by V. 

Let us denote by Slj C K.^, for j = 0, 1, 2, the set of resonant vectors of order exactly j, i.e., 

VLj := e such that = 2 - j}, with A^ {k e 1? such that ^ • fc = 0} = n Z^ 

Note that the sets VLj form a partition of M^, and that we have 
• f^o = {0}; 



^ e ill if and only if the geodesic issued from any x £ T in the direction ^ is periodic 



^ G if and only if the geodesic issued from any a; G in the direction ^ is dense in 



For each K gV such that rk(A) = 1, we define 

A-L e such that ^ • A: = for aU A: € A} , 

Ta := (A) /27rA, 
Ta-l := A-L/(27rZ2nA-L). 

Note that Ta and are two submanifolds of diffeomorphic to one-dimensional tori. Their 
cotangent bundles admit the global trivialisations r*TA = Ta x (A) and T*Ti^± = Ta-l x A-'-. 

For a function / on T^ with Fourier coefRcients {f{k))k£Z^ , and A e P, we shall say that / has 
only Fourier modes in A if f{k) = for k ^ A. This means that / is constant in the direction A^, 
or, equivalently, that a ■ d^f = for all a E A-^. We denote by i^(T^) the subspace of L^'(T^) 
consisting of functions having only Fourier modes in A. For a function / e i^(T^) (resp. a symbol 
a € S'°(T*T^)), we denote by (/)^ its orthogonal projection on L^(T^), i.e. the average of / along 
A^: 

(/)a E r^sp- («) A ^) E ^) ) ■ 

keA V fceA / 

If rk(A) = 1 and u is a vector in A-'- \ {0}, we also have 

1 



(/)a(x)= lim - / f{x + tv)dt. (7.1) 

T-^co I Jq 

In particular, note that (/)^ (resp. (a)^) is nonncgative if / (resp. a) is, and that (/)^ e "^""(T^) 
(resp. (a) J, S S°{T*T^)) if / G ^""(T^) (resp. a G S°iT*T^)). 

Finally, given / G L^(T^), we denote by nif the bounded operator on L^(T^), consisting in the 
multiplication by /. 
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7.2 Decomposition of invariant measures 

We denote by Al+(r*T^) the set of finite, nonnegative measures on T*T'^. With the definitions 
above, we have the following decomposition Lemmata, proved in [MaclO] or [AMll, Section 2]. 
These properties are given for general measures /x € M^(T*T^). Of course, they apply in particular 
to the measure fj. defined by Proposition 5.1. 

Lemma 7.1. Let fi G M^{T*T'^). Then ji decomposes as a sum of nonnegative measures 

M = /"|t2x{0} + M|t2x02 + X/ Ai|T2x(A^\{0}) (7-2) 

Ae-p,rk(A) = l 

Given /i e Ai^{T*T'^), wc define its Fourier coefficients by the complex measures on M^: 

jx{k, ■) := I — pL{dx, •), k IE Z. 

One has, in the sense of distributions, the following Fourier inversion formula: 

Lemma 7.2. Let /i e A^+(T*T^) and KeV. Then, the distribution 

^ ^ik-x 

keA 

is in M+{T*T'^) and satisfies, for all a e 'S'^{T*T'^), 

((a')a ' '^)a^(T'T2)x^'(t*t2) ~ {f-T {'^) id M(T'i:^),'g^(T-'i^) ■ 

Lemma 7.3. Let /i g M~^{T*'f'^) he an invariant measure. Then, for all A E V, /^|t2x(a^\{o}) 
also a nonnegative invariant measure and 

A'lT2x(Ai\{0}) = (m)a lT2x(Ai\{0})- 

Let us now come back to the measure fi given by Proposition 5.1, which satisfies all properties 
listed in Proposition 6.1. In particular, this measure vanishes on the non-empty open subset of 
given by {b > 0} (see Item 4 in Proposition 6.1). As a consequence of Proposition 6.1, and of the 
three lemmata above, this yields the following lemma. 

Lemma 7.4. We have ^ = EAeP,rk(A)=i MlT2x(Ai\{o})- 

As a consequence of Proposition 6.1, wc have indeed that the measure fi is supported in {|^| = 1}, 
which implies /i|T2x{o} = 0. In addition. Lemma 7.3 applied with A = {0} implies that ^|t2x02 
constant in a; - and thus vanishes everywhere since it vanishes on {b > 0}. 

Remark 7.5. Since the measure ^ is supported in {|^| = 1} (Proposition 6.1, Item 1), we have 

MIt^xA^ = M|t2x(A^\{0}) 

(which simplifies the notation). 

As a consequence of these lemmata and the last remark, the study of the measure /i is now 
reduced to that of all nonnegative invariant measures /iIt^xA^ with rk(A) = 1. 

The aim of the next sections is to prove that the measure /x|t2xA^ vanishes identically, for each 
periodic direction A^. 
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7.3 Geometry of the subtori Ta and T^x 

To study the measure /i|T2x(AX\{o})i need to describe briefly the geometry of the subtori Ta and 
Tjy± of T^, and introduce adapted coordinates. 
We define xa the linear isomorphism 

XA : A-L X (A) ^ R2 : {s,y)^s + y, 

and denote by xa • r*A-'- x T* (A) r*R^ its extension to the cotangent bundle. This application 
can be defined as follows: for {s,(t) e r*A-'- = A-'- x (A-'-)* and (y,?^) £ T* (A) = (A) x (A)*, we can 
extend cr to a covector of vanishing on (A) and 1] to a covector of vanishing on A-'- . Remember 
that we identify (R^)* with R^ through the usual inner product; thus we can also see a as an element 
of A-'- and as an element of (A) . Then, we have 

Xa(s, a, y, 7?) = (s + y, a + 7?) e T*R^ = R^ x (R^)* . 

Conversely, any ^ e (R^)* can be decomposed into f = cr + /y where a g A-^ and t] E (A). We 
denote by Pa the orthogonal projection of R^ onto (A), i.e. = -q. 

Next, the map xa goes to the quotient, giving a smooth Riemannian covering of by 

: Ta^ X Ta ^ : (s, y) ^ s + y. 

We shall denote by tta its extension to cotangent bundles: 

*A : r*TA-L X r*TA ^ t*t2. 



As the map tta is not an injection (because the torus Ta-l x Ta contains several copies of T^), we 
■oducc its degree pA, w 
Then, the application 



introduce its degree pA, which is also equal to ^"'^"^^'Ij^'^'* 



Tau := o XA, 

v/PA 



defines a linear isomorphism Lj^Q^(R^) Lif^^{K^ x (A)). Note that because of the factor -^^i 

Ta maps i^(T^) isometrically into a subspace of L'^{Tj^± x Ta). Moreover, Ta maps L^(T^) into 
^^(Ta) C L^(TAi X Ta), since the nonvanishing Fourier modes of m € L^(T^) correspond only to 
frequencies fc e A. This reads 

TA«(s,y) = -^u{y) for (.s,y) e Ta^ x Ta. (7.3) 

\/PA 



Since xa is linear, we have, for any a S '^°°(T*R^) 

TaOpJa) -Op„(aoxA)TA, (7.4) 

where on the left Op^ is the Weyl quantization on (A.l), and on the right Op^ is the Weyl 
quantization on A^ x (A). Next, we denote by Op^ and Op^ the Weyl quantization operators 
defined on smooth test functions on T*A-'- x T* (A) and acting only on the variables in T*A-'- and 
T* (A) respectively, leaving the other frozen. For any a e '^^°°(T*A-^ x T* (A)), we have : 

Op„(a) = Opf o Op^(a) = Op;^ o Opf (a). (7.5) 

Now, if the symbol a £ '^#'^°°(T*T^) has only Fourier modes in A, we remark, in view of (7.3), that 
ao TTA does not depend on s G Ta^ . Therefore, we sometimes write ao iTA^a, y, rj) for ao tta(s, cr, y, rj) 
and (7.4)- (7.5) give 

Ta Op,,(a) = Op;^ o Op^ (a ° ^a)Ta = Op;^(a o n^ikD,, •, •))Ta. (7.6) 

Note that for every a € A-'-, the operator Op^ (a o tta (cr, •, •)) maps T^(Ta) into itself. More precisely, 
it maps the subspace Ta(Ta('^^)) i^ito itself. 



21 



8 Change of quasimode and construction of an invariant cut- 
off function 



In this section, we first construct from tlie quasimode Vh another quasimode Wh, that will be easier 
to handle when studying the measure fi\j'2y,/^±. Indeed Wh is basically a microlocalization of Vh in 
the direction at a precise concentration rate. 

Moreover, we introduce a cutoff function Xhi^) — XhiUi^); well-adapted to the damping coeffi- 
cient b and to the invariance of the measure ^|t2xA^ in the direction (this cutoff function plays 
the role of the function xib/h) used in [BH07] in the case where b is itself invariant in the direction 
A^). Its construction is a key point in the proof of Theorem 2.6. 

Let X G 't^^iM.) be a nonnegative function such that x = 1 in a neighbourhood of the origin. We 
first define 

Wh:= Op,, (^X (^^))^'" 

which implicitely depends on a G (0, 1). The following lemma implies that, for 6 and a sufficiently 
small, Wh is as well a o(/i^+'')-quasimode for PI^. 

Lemma 8.1. For any a > such that 

S+^ + a<^, 3a + 2(5<l, (8.1) 

we have 

|lnV|U.(T2)=o(/i2+^). 

As a consequence of this lemma, the semiclassical measures associated to Wh satisfy in particular 
the conclusions of Proposition 6.1. Moreover, the following proposition implies that the sequence Wh 
contains all the information in the direction A^. 

Proposition 8.2. For any a E '^^{T*T'^) and any a £ (0,3/4) satisfying the assumptions of 
Lemma 8.1, we have 

Next, we state the desired properties of the cutoff function Xh- The proof of its existence is a 
crucial point in the proof of Theorem 2.6. 

Proposition 8.3. For 5 — 8e, and e < -g, there exists a satisfying (8.1), such that for any constant 
Co > 0, there exists a cutoff function Xh € ^°°(T^) valued in [0, 1], such that 

1. Xh ~ Xhiy) does not depend on the variable s (i.e. Xh h.^ -invariant) , 

2. \\{\~xt)wh\\L^iT-)^o{\), 

3. b< coh on supp(x^), 

4- WdyXh^hh^T^) = o(l), 

5- Wd^XhWhWmT^) - o(l). 

Note that the function Xh implicitely depends on the constant cq, that will be taken arbitrarily 
small in Section 10. 

In the particular case where the damping function b is invariant in one direction, this proposition 
is not needed. In this case, one can take as in [BII07] xt ~ d-dimensional torus, this 

cutoff functions works as well if b is invariant in d — 1 directions, and an analogue of Theorem 2.6 
can be stated in this setting. Unfortunately, our construction of the function xt (see the proof of 
Proposition 8.3 in Section 13) strongly relies on the fact that all trapped directions are periodic, and 
fails in higher dimensions. 

We give here a proof of Lemma 8.1. Because of their technicality, we postpone the proofs of 
Propositions 8.2 and 8.3 to Sections 12 and 13 respectively. 
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Proof of Lemma 8.1. First, we develop 



Pr 



W - Pi Op. (x (^)) V, ^ Op, (x (^) ) Pivn + ^h 



Vh, 



since Pq and Op, (^^^)) ^'"^ both Fourier mutipliers. We know that 



Op, X 



< 



L2(T2) 



\\PS-,WL-HT^)=o{h'+'). 



It only remains to study the operator 



= Z/.1-0P, ( dyb x' (^-^))+OciLMh'''-''^) 



^Op, (^x(^^^^^j ... ^^nyv-^ \ h-' J 
according to the symbolic calculus. Moreover, using Assumption (2.12), we have 



< Ch 



l-e 



The sharp Carding inequality applied to the nonnegative symbol 

2 

then yields 



^2^2(l-e) _ 



Op, C^6^(^-) 



and hence 



Op. ( X' I 



L2(T2) 



< C2(52(l-')z;,,z;,)i.(T2) + 0(/il-"). 



When using the inequality J /-'^ "di^ < (J fdv^^ ^ for nonnegative functions (with dv = \vh{x)\'^dx), 
we obtain 

Combining this estimate together with (8.3) and (8.4) gives 



(8.3) 



(8.4) 



ih 



b, Op, X 



Vh 



L2(T2) 



Coming back to the expression of Pj^Wh given in (8.2), this concludes the proof of Lemma 8.1. □ 



9 Second microlocalization on a resonant afRne subspace 

We want to analyse precisely the structure of the restriction A'It2x(A-|-\{o})j using the full information 
contained in o(/i^^'^ )-quasimodes like w, and w,. 

From now on, we want to take advantage of the family w, of o(/i^''"'')-quasimodcs constructed in 
Section 8, which are microlocaliscd in the direction A^. Hence, we define the Wigner distribution 
y^rh g ^'(r*T2) associated to the functions Wh and the scale h, by 

a)(50),,so = (Op,(a)u;,, «;,)^.(T2) for all a £ S\T*T^). 
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According to Proposition 8.2, we recover in the limit /i — >■ 0, 

for any a G '^;?°{T*T'^) (and a satisfying (8.1)). 

To provide a precise study of ii\j2~^f^±, we shall introduce as in [MaclO, AMll] two-niicrolocal 
semiclassical measures, describing at a finer level the concentration of the sequence Vh on the resonant 
subspace 

= {^e such that P^S, = 0}. 

These objects have been introduced in the local Euclidean case by Nicr [Nie96] and Fermanian- 
Kammerer [FKOOb, FKOOa]. A specific concentration scale may also be chosen in the in the two- 
microlocal variable, giving rise to the two-scales semiclassical measures studied by Miller [Mil96, 
Mil97] and Fermanian-Kammerer and Gerard [FKG02]. 

We first have to describe the adapted symbol class (inspired by [FKOOa] and used in [AMll]). 
According to Lemma 7.3 (see also Remark 7.5), it suffices to test the measure ^\j2xA^ with functions 
constant in the direction A^ (or equivalently, having only x-Fourier modes in A, in the sense of the 
following definition). 

Definition 9.1. Given A S "P, we shall say that a € ii a = a{x,^,ri) e '^°°{T*T'^ x (A)) and 

1. there exists a compact set Ka C T*T'^ such that, for all rj E (A), the function (x, ^) h-> a{x, ^, ?/) 
is compactly supported in Ka', 

2. a is homogeneous of order zero at infinity in the variable rj e (A); i.e., if we denote by Sa := 

n (A) the unit sphere in (A) , there exists i?o > (depending on a) and Ohom G ^c°° {T*T'^ x §a) 
such that 

a{x, V) = flhom (^x, , for \rj\ > Ro and (.t, ^ G T*T^; 

for ?7 7^ 0, we will also use the notation a{x, ^, ooi]) := flhom (^x, ^, 

3. a has only x-Fourier modes in A, i.e. 

a{x,£„v) ^ V ^—a{k, £,,?]). 

^ — ' ZTT 

feeA 

Note that this last assumption is equivalent to saying that a ■ dxa = for any a E A-^. We denote 
by S\ the topological dual space of S*^. 

Let X S "^(""(li) be a nonnegative function such that x = 1 hi^ a neighbourhood of the origin. Let 
R > 0. The previous remark allows us to define, for a E S\ the two following elements of Sj^' : 

.= (H.-Vfl-.(^))<.(.,«.i|i)\ , (9.1) 

(^A^ := {W\X (^) a L e, ^) ) . (9.2) 



@'(T*T2),<g'_?=(T"-T2) 

In particular, for any R > and a E S\, we have 

/w'\a L e, ^) \ = (l^r . a) , + {Wl,, a) . (9.3) 

The following two propositions are the analogues of [FKOOa] in our context. They state the 
existence of the two-microlocal semiclassical measures, as the limit objects of W'^' and ^. 
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Proposition 9.2. There exists a subsequence {h,Wh) and a nonnegative measure S A^ + (T*T^ x 
§a) such that, for all a G S\, we have 

lim \im/w^'^,a\ ^ = /i^^,ahom [x,(,, y- . . 

R^ooh-^0\ ■ ' Sj^ \ \ I'?! / / 7M(T*T2xSa),'^0(T*T2xSa) 

To define the limit of tlie distributions Wj^ j^, we need first to introduce operator spaces and 
operator- valued measures, following [Ger91]. Given a Hilbert space H (in the following, we shall 
use H = L^(Ta)), we denote respectively by IC{H), C^{H) the spaces of compact and trace class 
operators on H . We recall that they are both two-sided ideals of the ring C{H) of bounded operators 
on H. We refer for instance to [RS80, Chapter VI. 6] for a description of the space C^{H) and its 
basic properties. Given a Polish space T (in the following, we shall use T = T*Ti^±), wc denote by 
A^ + (T; C^{H)) the space of nonnegative measures on T, taking values in C^{H). More precisely, we 
have p e A^+(T; £^(iJ)) if p is a bounded linear form on '^^(T) such that, for every nonnegative 
function a e 'if^{T), {p 

' ^)m{t) V{t} ^ ^^i^) ^ nonnegative hermitian operator. As a consequence 
of [RS80, Theorem VI. 26], these measures can be identified in a natural way to nonnegative linear 
functionals on '^^(T; /C(iJ)). 

Proposition 9.3. There exists a subsequence {h,Wh) and a nonnegative measure 

PAeM+{T*T^^;C\L^TA))), 
such that, for all K G 'i^^{T*Ta^;JC{L'^{Ta))), 



\\m{K{s,hDs)TAWh,TA_Wh)i^2n . ■L'^n^))=^T: { / K{s,a)pK{ds,da)) . (9.4) 

/l— >-0 ^ A ' ^ " I Irp^j.^ 

Moreover (for the same subsequence) , for all a G S\, we have 

hm \im{W^^^,a) , =tr\ [ Op^ (a(^A(a, y, 0), 7,))pA(ds, da) I . (9.5) 

In the left hand-side of (9.4), the inner product actually means 
{K{s,hDs)TAWh,TAWh)L2^^r^^L2(^j^^') 



e' 



In the expression (9.5), remark that for each a e A-'-, the operator Op^ (a(7rA(cr, y, 0), 77) is in 
£(L^(Ta)). Hence, its product with the operator pA{ds,da) defines a trace-class operator. 

Before proving Propositions 9.2 and 9.3, we explain how to reconstruct the measure A^It^xA^ 
from the two- microfocal measures and pA- This reduces the study of the measure p to that of all 
two-microlocal measures and pA, for A G V. 

We denote by A^+(T) the set of compactly supported measures on T, and by (•, •)^ c^oi^t) 
the associated duality bracket. 

Proposition 9.4. For all a e '^^^ {T*T'^) having only x-Fourier modes in A (i.e. for all a G ^a 

independent of the third variable rj £ (A) ), we have 



(a*: '3^)a^(T*T2)X0(T*T2) — ('^^J '3^)7V1(T*T2xSa),'^^'c°(^*T^xSa) 1 / "^ao-S-A (''')PA (^s, dcr) > , (9.6) 

Jt*t.j_ I 
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and 

(mIt^xA^i 0)7^1 (T'T2),'i!fi' (T*T2) ^ ('^'^lT2xA^xSA'")7H(T*T2xSA)Xe°(T*T2xSA) 

+ tr<; / mQo#A(cr)pA(ds, da) > , (9.7) 
{Jt'T^^ J 

where for a G A^, TOaoirA(o') denotes the multiplication in L^(Ta) by the function y 1— > a o TT/^{a,y). 

Moreover, we have G A^+(T*T^ x Sa) a?fi Pa G (T'*irAi ; ^^(Ta)) ('i.e. &ot/i measures are 
compactly supported) . 

Formula (9.7) follows immediately from (9.6) by restriction. By the definition of the measure pA, 
we see that it is already supported on x A-*- (see expression (9.2)). 

The end of this section is devoted to the proofs of the three propositions, inspired by [FKOOa, 
AMll]. 

Proof of Proposition 9.2. The Calderon-Vaillancourt theorem implies that the operators 

°- (0 - (^)) « (^'^^ ^)) (0 - (¥)) 

are uniformly bounded as /t, — > and R — > +cx). It follows that the family W'^^ is bounded in 5a', 
and thus there exists a subsequence (/i, wu) and a distribution jl^ such that 



lim lim ( T^p''^, a) = (/i^, a (x, ^, r/)) 1 , 1 . 



Because of the support properties of the function x, we notice that (/i^, a) „i , „i = as soon as the 

^ ' -^A '-^A 

support of a is compact in the variable 77. Hence, there exists a distribution G Si'{T*T'^ x Sa) 
such that 



V]\J I @'(T*T2xSa),'i?~(T*T2xSa) 

Next, suppose that a > (and that ^\ — x is smooth). Then, using [AMll, Corollary 35], and 
setting 



c 

< —. 



£(L2(T2)) 



^ Rh J J V 

there exists C > such that for all h < ho and i? > 1, we have 

As a consequence, we have, 

/W^'^,a\ > II Op;,(6^)u;,,||^2(T2) - plk/i|li2(T2), 

^ ' -^A '■-'a ^ 

SO that the limit ( J^"^,ahom -rr ) / is nonnegative. The distribution i'^ 

\ \ ' / / ®'(T*T2xSa),'^~(T*T2xSa) 

is nonnegative, and is hence a measure. This concludes the proof of Proposition 9.2. □ 

Proof of Proposition 9.3. First, the proof of the existence of a subsequence (/i, Wh) and the measure 
Pa satisfying (9.4) is the analogue of Proposition 5.1 in the context of operator valued measures, 
viewing the sequence Wh as a bounded sequence of L^(TAi ; ^^(Ta)). It follows the lines of this 
result, after the adaptation of the symbolic calculus to operator valued symbols (or more precisely, 
of [Ger91] in the semiclassical setting). 
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Second, using the definition (9.2) together with (7.6), we have 



a he, 4, —r- Wh,Wh 



Rh V /i 



L2(T^^xTa) 



Hence, setting 



we obtain 

<W^fl,A,a)si'.5i = (Opf oOp^ (a5j,A(^,2/,^/))TA^«/^,TAi(;„ 



L2(T^_lxTa) 



We also notice that Op^ (a^ a) ^ /C(£^(Ta)), for any cr G A-"- since a '^^s compact support with 
respect to 77. Moreover, for any R > fixed and a E S\, the Calderon-Vaillancourt theorem yields 



Opf (<,a) = Op;^ (a'A,A) + hB 

for some B G £(L^(Ta)), uniformly bounded with respect to h. Using (9.4), this implies that for 
any R > fixed and a G ^a, we have 



hm (Ty^,A, SI - tr M Op^ «A)PA(ds, da) . 

Moreover, we have 

hin Op^ {a%^) = Op^ (alL.A) = Op^ (a(^A(a, y, 0), r;)) , 

in the strong topology of '^^{T*Tf^±; C(L'^{Ta)))- This proves (9.5) and concludes the proof of 
Proposition 9.3. □ 

Proof of Proposition 9.4- Taking a G Sj^, independent of the third variable rj G (A) gives 

(W , a (a;, ^))^,^j,,^2)^.^oo(y.T2) ~^ {f^\T^xA-^TO,)M{T'T'^),V°{T*T^) ^ 

together with 

■'a '■'a 

(according to Proposition 9.2) and 

(Wra^o) I, , ^ tT \ / Op^ {a{TTAi(T,y,0)))pA{ds,da)\ =tT < / maoirA{'^)PA{ds,da) } , 

(according to Proposition 9.3). Now, using the last three equations together with Equation (9.3) 
directly gives (9.7). 

As both terms in the right hand-side of (9.7) are nonnegative measures and the left-hand side 
is a compactly supported nonnegative measure, this implies that i/^ and pA are both compactly 
supported. □ 
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10 Propagation laws for the two-microlocal measures and 

PK 

In this section, wc study the propagation properties of and piy. The key point here is the use of 
the cutoff function introduced in Proposition 8.3. 

We wiU use repeatedly the following fact, which follows from Item 2 in Proposition 8.3: if A is a 
bounded operator on i^(T^), we have 

{Awh,Wh)L'^(T'^) = {AXhWh, XhWh) ^2fj2) + ll^ll£(L2) o(l). (10.1) 

To simplify the notation, we shall write Ac^ji for Xh-^Xh- 
10.1 Propagation of 

We define for (cc, 77) G T*T'^ x (A) and t e M the flows 
generated by the vector field ^ • dx and, for rj ^ 0, 

generated by the vector field ■ dx- With these definitions, we have the following propagation laws 
for the two-microlocal measure v^. 

Proposition 10.1. The measure is (j)^- and (jy^.- invariant, i.e. 

{f^).,v^ = v^ and ^v^, for every t eM.. 

The key result here is the additional "transverse propagation law" given by the flow The 
measure not only propagates along the geodesic flow 0°, but also along directions transverse to 

Proof. Fix a ^ S\. The computation done in (6.2) is still valid replacing a by ^1 — x { ^^R}f ^ ) ^ "IT 

since it only uses the fact that Op/j ((l ^ X i ^^Rh )) ^ ^' "^T^) ) bounded and that HP^'^w^ 11^2(^2) = 
o{h) and {hwh,Wh) ]^2(j2) = o{l). This yields 



\im(w^-^,^-dxa] 



lim(H"-,f-&{(l-i 



Rh V /i 



®'(T*T2)Xc°°(T*T2) 



and hence, in the limit R — > +cx), we obtain 

\''1\J / A1(T*T2xSA),<g'^"(T*T2xSA) 

Replacing ahom by Ohom ° ^'^d integrating with respect to the parameter r gives (0^)*i^^ = i^^, 
which concludes the first part of the proof. 

Second, to prove the 0]--invariance of we compute 
/i/^,^-9,ahom = hm \iTnlw^^,^^-dxcS . (10.2) 

\ I'^l V 1^/1/ / A^(T'T2xSa)X,"(T*T2xSa) R^^h^O\ |7/| I S]^' ,Sl 
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Setting 



and 



we have the relation 



• h J 



h.A V 



R ' 



v\ 



■ dxtt 



where Aa = dy is the laplacian in the direction A. 



(10.3) 



Lemma 10.2. For any given cq > and R > 0, we have 

{[Aa, A^']wh,Wh)mT^) = {[^A, A^o,h\wh,Wh)L^(^T^) +o(l). 

We postpone the proof of Lemma 10.2 and first indicate how it allows to prove Proposition 10.1. 
We now know that 



A J]_ p. I c: 

\m 



lim \im--{[AA,A^^ f^]wh,Wh)L^T^)- 



Recall that a G S\ implies that a has only x-Fourier modes in A, i.e. Pa£, • dxCL = ^ ■ d^a. We have 
also assumed in this section that h has only x-Fourier modes in A. As a consequence, we have 



'-([Aa, Af^jJw/,,u;/0L2(T2) = --([A, Af^ «7h)i2(T2) 



2h? 



L2(T2) 



(10.4) 



Developing the last expression of (10.4), we obtain 
i 

2h^ 



(10.5) 



Since ^ is bounded in £(i^(T^)), its adjoint is also bounded so that the first two terms in 

the last expression vanish in the limit /i — > 0, using 11^2(^2) = o(/i^). To estimate the last two 

terms, we use again the boundedness of A^ and (A^)* and write 

I «,/^^^,Ml2(t2) I ^ ll^""!! IIx^^«/..||l2(t2) < 2c„h p^ii, 

according to Item 3 in Proposition 8.3. It follows that 
1 , . ^ 1 



lim sup 



2h 



Coming back to the expression (10.2), we obtain 



V 



Vl/ / X(T*T2xSa),'^°(T*T2xSa) 



< 2coSupP^ | 



< 2co supP^ 



and since Cq was arbitrary, 



\m 



1^1/ / >l(T*T2xSA),'if°(T*T2xSA) 



= 



Replacing ahom by Ohom o 4>l ^^'^ integrating with respect to the parameter r gives ((^^)*i^^ = ly^, 
which concludes the proof of Proposition 10.1. □ 
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Proof of Lemma 10.2. We are going to show that 

([Aa, w/0l2(t2) = {[Af„A^]coMWh,Wh)mT^) +o(l). (10.6) 

Then, using the fact that [Aa, A^'] is a bounded operator (its symbol is — x (t?')) ]ri\''^^^^^^ ^' ^))' 
together with (10.1), this is also {[Aa, A-'^]wh,Wh)L^(T2-f +o(l). 

To prove (10.6), we develop the difference [Aa, ,J — [Aa,A^]co,/i as 

[Aa,<,,J - [Aa,A%„m = Kxt] + xi^A« [dlx^i] . (10.7) 

Then, writing 
we have 

{[dy,Xh] A^Xh'Wh,Wh)^2^.^2-j = {A^Xh Wh, dlxt Wh) ^2^j2^ + {dy o A^ XhWh,'2.dyxtwh) ^^(j^) ■ 
Recalling that the operator dyoA^ is bounded, and using Items 4 and 5 in Proposition 8.3, we obtain 

{[dlxh] A\twH,WH) ^2^^,^ < C\\d'yXtwh\\LHT^-)++C\\dyX^iw,]\L2^j2^^o{l). 

The last term in (10.7) is handled similarly. This finally implies (10.6) and concludes the proof of 
Lemma 10.2. □ 

10.2 Propagation of 

We denote by {lj-'^, e'J^^)jgN the eigenvalues and associated eigcnfunctions of the operator — Aa = ~dy 
forming a Hilbert basis of L^(Ta). We shall use the projector onto low frequencies of — Aa, i.e., for 

any co e M_|_, the operator 

which has finite rank. 

We have the following propagation laws for the two-microlocal measure pA- 

Proposition 10.3. 1. For any K <E {T*Ta±;IC{L'^{Ta))) , independent of s (i.e. K{s,a) = 
K{a)), and any lo > 0, we have 

tr J / [Aa, U%K{a)U%] pAids, da) I = 0. 

2. Moreover, defining 

A/a:- / pA{ds,d<j)eC\L'-{TA)), 

we have 

[Aa,Ma] =0. 

Remark that for any a G A^, the operator 

[AA,nx/^(a)nx] = nx[AA,/v(a)]nx, 

has finite rank, so the right hand-side of Item 1 is well-defined. Note that the definition of Ma has 
a signification since pA has a compact support, according to Proposition 9.4. 

The commutation relations of Items 1 and 2 in this proposition correspond to propagation laws 
at the operator level. They are formulated here in a "derivated form" , which, for Item 2 for instance, 
is equivalent to 



e^'^'^'^MAe-''"^'' = Ma, for aU t e 



in the "integrated form" . 



30 



Proof. For K G 'g;°°(A^; /C(L2(Ta))) (in other words K € t^^{T*Tj^±;lC{L'^{TA))) independent of 
s £ T^i ) , we denote 

and wc note that K'^ is also in '^^^{K^\JC{L'^{Tfs.)))- Hence, we have 



tr 



/ [AA,ni(X(a)nX] PA[ds,da) \ = - lim ([-Aa, A'"(/iZ?,)]rAu;,„ Taw^O^.^ 



)) 



To show that this limit vanishes, wc proceed as in lines (10.4), (10.5) and in the subsequent calcula- 
tion, replacing the operator by K'^{hDs). 

With the notation Aa = dy and Aa^ = d^, we first note that 

{[-Aa, K'^ {hDs)]TAWh,TAWh) L2^T^^.L^r^-^-^ = i[-A,K'^{hD,)]TAWh,TAWh)L2(^j^^.L2(^j^)) , 

since A Aa + Aa^ and since [A/^±,K'^{hDs)] = 0. As a matter of fact, K'^{hDs) = Opf^{K'^{a)) 
and Aa-l = —h~'^ Op^(|(Tp) are both Fourier multipliers. 

The following lemma is proved the same way as Lemma 10.2 

Lemma 10.4. For any given cq > 0, we have 

{[Aa, K'^ icr)]TAWh,TAWh)L^{T^) = {[Aa, K'^^ji{hDs)]TAWh,TAWh)L^{T2) +0(1). 

Here K^^j^QiD,,) means Xh^"^ {hD,)xt. 
Writing 

-h'^A = TaPI'TI - ihb o tta, 

we have 

{[-A,K%,h{hDs)]TAWh,TAWh) ^^^^^^.^^^^^^^^ 

= ^ {K,,h{hD,)TAWh,TAPtwh)^^^^^^.^^^^^^^ - ^ {K'^^,^fXhDs)TAPtwh,TAWh)^,^^^^^.^^^^^^^ 
+ 1 {Kuh{hDs)TAWh,TA{hwh))^.,^^^^,^,^^^^^ + ^ {K^^j,{hDs)TA{hwh),TAWh)^,^^^^,^.,^^^^^ . 
It follows, as in (10.5), that 

lim^upl {[-A,K'^^^ ,^{hDs)]TAWh,TAWh)^^^^^^.^^^^^^^ \ < 2co\[K\[ 

and since co was arbitrary, we can conclude that 

\im{[AA,K^{<j)]TAWh,TAWh)mT^) = 0, 

which concludes the proof of Item 1. 

Item 1 gives, for all K E IC{L'^{Ta)) constant (which is possible since pA{ds,da) has compact 
support). 



= tr J / [Aa, K'^]pA{ds, da) \ = tr J [Aa, A'"] / PA{ds, da) \ = tr {[Aa, A'"]A/a} • 

Using that tr (AB) = tr {BA) for all A £ £^ and B E C together with the linearity of the trace 
(see [RS80, Theorem VI.25]), we now obtain, for all K G IC{L'^{Ta)), and aU w > 0, 

= tr {[Aa, nXA'HX]il/A} - tr {if HX[Aa, Ma]HX} . 

Consequently, we have for all a; > 0, HX[Aa, Ma]HX = (see [RS80, Theorem VI.26]). Letting w go 
to +00, this yields [Aa, Ma] = and concludes the proof of Item 2. □ 



31 



11 The measures and vanish identically. End of the 
proof of Theorem 2.6 

In this section, wc prove that both measures and pA vanish when paired with the function {b) j^. 
Then, we deduce that these two measures vanish identically. In turn, this implies that /i|T2xA-L = 0, 
and finally that ^ = 0, which will conclude the proof of Theorem 2.6. 

Proposition 11.1. We have 

(j^^^It^xA^xSa' Wa)a^,(t*t2xSa),"^''(t*t2xSa) "''^'^ tr{m(h)^MA} = 0. 

As a consequence, we prove that p\ and i'^It^xA^xSa vanish. 
Proposition 11.2. We have pA = and i^^It^xA^xSa ~ ^- Hence /i|T2xA-i- = 0. 

This allows to conclude the proof of Theorem 2.6. Indeed, as a consequence of the decomposition 
formula of Proposition 9.4, we obtain, for all A £ V, such that rk(A) = 1, /x|t2xa-|- = 0. Using the 
decomposition of the measure given in Lemma 7.1 together with Lemma 7.4, this yields fx = 
on T^. This is in contradiction with /i(T*T^) = 1 (Proposition 6.1), and this contradiction proves 
Theorem 2.6. 

We now prove Propositions 11.1 and 11.2 
Proof of Proposition 11.1. First, (4.22) implies that {bvh,Vh)L'^{T^) — >■ 0, and hence 

(M> ^)a^4T'T2)X°(T*T2) = 0- 

Then the decomposition given in Lemma 7.1 into a sum of nonnegative measures yields that, for all 
A e 7^, 

(mIt^xA-Lj ^)^^(T*T2),<g'0(y.T2) ~ 0, (11-1) 

since h is also nonnegative. Lemmata 7.2, 7.3 and 7.4 (see also Remark 7.5), then give 
(Ml 

= {^A^^xA^^^ Ma{T'T'^),'i^0(T'V) (11-2) 

where the function (6)^ is also nonnegative. The decomposition formula of Proposition 9.4 into the 
two-microlocal semiclassical measures then yields 

(A^It^xA^j (^)a)7Wc(T*T2),'^0(T*T2) = ('^'^It^xA^xSa' a) Ma(T^T^x%^_),V°{T'T^y^^f,) 

+ tr < / mi^i,) pA{ds,da) \ . 

Besides, the measure j/^|t2xA-LxSa is nonnegative, hence (i^^It^xA^xSa' (^)a)>i^(t*t2xSa)."»'0(t-t2xSa) - 
0. Similarly, p\ g 7M + (T*Ta-l ; £^(Ta)) and the operator m(^h) ^ € '^(^^(Ta)) is selfadjoint and non- 
negative, which gives tr | J^,,^ ^ m(j^^ ^pjy{ds.,da)^ > 0. Using (11.1) and (11.2), this yields 

('^^It^xA^xSai (^)a)a1c(T*T2xSa),'^0(T*T2xSa) ^ ^' 

and 



W m(^i,^^pA{ds,da) S = 0. 



In this expression, the operator rnij^'^^ does not depend on (s, ct), so that 



= tr < TO(b)^ / pA{ds,da) ^ = tr{m(6)^ A/a}, 
which concludes the proof of Proposition 11.1. □ 
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Proof of Proposition 11.2. Let us first prove that pA = 0. We reeah that the operator Ma is a 
selfadjoint nonnegative trace-class operator. Moreover, Proposition 10.3 implies that the operators 
Ma and Aa commute. As a consequence, there exists a Hilbert basis {e'/^jtz^ of L^(Ta) in which 
Ma and Aa are simultaneously diagonal, i.e. such that 

-^Ae^A = ^A^A^ a,nd Mac^ = 7a ^A' 

where (7A)iGN are the associated eigenvalues of Ma. In particular, we have 7^^ > for all j S N (and 
7^ G £^). Note that the basis (ey^)jgN is not necessarily the same as the basis (ey^)jgN introduced in 
Section 10.2. 

Using Proposition 11.1, together with the definition of the trace (see for instance [RS80, Theorem 
VI. 18]) we have 

= tr{m(b)^MA} = J2 {^{b)^MAii,ei) = E^A (WAel,el) . 

Since all terms in this sum arc nonnegative (because both 7^ and (&)^ are), we deduce that for all 

Suppose that 7^ 7^ for some j E N. Then, I (6)^ e^^, e^^ ) =0 where (fe)^ is nonnegative and 

V / L^{Ta) 

not identically zero on Ta. This yields e^^ = on the nonempty open set {{b)A > 0}. Using a unique 
continuation property for eigenfunctions of the Laplace operator on Ta , we finally obtain that the 
eigenfunction e^^ vanishes identically on Ta- This is absurd, and thus we must have 7^ = for all 
j e N, so that Ma = 0. Since pA € A^+(T*Ta-l;£^(Ta)), this directly gives pA = 0. 

Next, we prove that J/^ = 0. This is a consequence of the additional propagation law of i^^ with 
respect to the flow (j)]. (see Section 10.1). Indeed the torus Ta has dimension one, = 
(according to Proposition 10.1) and, using Proposition 11.1, vanishes on the (nonempty) set 
{(6)a > 0} X R2 X §a (with > 0} clearly satisfying GCC on Ta). Hence, = 0. 

To conclude the proof of Proposition 11.2, it only remains to use the decomposition formula (9.7) 
which directly yields /iljaxA^ = 0. □ 

12 Proof of Proposition 8.2 

In this section, we prove Proposition 8.2. For this, we consider two-microlocal semiclassical measures 
at the scale h°' . The setting is close to that of [FK05]. 

We shall see that the concentration rate of the sequence Vh towards the direction A-'- is of the 
form for aU a < 2±i. 

First, Lemma 7.3 yields p\t^xA^ — {(Aa It^xA^ (see also Remark 7.5), i.e. 

and it suffices to characterize the action of /ij^axA^ on A^-invariant symbols. Recall that, for all 

(I^t'A M{T'T^),'^°{T'T^) ~ 1™ (Op/i(a)^^/i, W/j)^2(T2') ■ 

In this section, the assumption Vb G 'Tf°°{T'^) is used in an essential way for the propagation 
result of Lemma 12.2 below. Like in (9.1) and (9.2), let us define : 



.:=(r*,(l-xf!#))a(..«,^|i)) . (12.1) 



®'(T*T2)X~(T*T2) 



@'(T*T2),<g'~(T*T2) 
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for a e S\. 

We take R = R{h) /i^^i"") for some a € (0, 1), so that Rh = /i". The proof of Proposition 9.2 
apphes verbatim and shows the existence of a subsequence {h,Vh) and a nonnegative measure £ 



hm ( , a ) , , = ( , flhom ( x, ^, ) ) 
/,_,o\ ^^(") /si', Si \ V HJ/j 



1^1/ / A^(T*T2xSA)Xe°(T*T2xSA) 

Proposition 12.1. Let R{h) = /i-^i""' mf/i a < T/ien 

'^a It2x(A-I-\{0})xSa = 

The proof of Proposition 12.1 relies on tlie following propagation result. 
Lemma 12.2. For a < the measure is (j)^- and (f>^ -invariant, i.e. 

{4P^),v^ = vt and {4>\),u^^u^, for every t eM.. 

The proof is very similar to that of Proposition 10.1 but docs not use Assumption (2.12). 

Proof. The proof of (/)"-invariance is strictly identical to what has been done for Proposition 10.1 and 
thus we focus on the 0^-invariance. Equation (10.5) still holds with R{h) ~ h~^^~"\ now reading 

\V\ /si 'si ^ /L2(T2) 2h^ V " /L2(T2) 



2h V /L2(T2) 2h V /L2(T2) 

where was defined in (10.3). Using \\Pl^Vh\\L^(T2-^ = o{h^^^) together with the boundedness of 
A^C*) , it follows that 



lim ( v!l:t\ i-p-i ■ dr 



Recall from (4.22) that ||\/6u;i||^2(-t2-) = o{h 2 ). In addition, it follows from standard microlocal 
calculus that 

[A^^'^\Vb]^OciL-){Rih)-'). 

We can thus write 

^, ^=o{l)-^{A^'i'^^Vbv,,Vbv,)^,^^^^ + l-{Vb[A^^^^^ 

\ I /I / Sa ,Sa 

+ ^([A^''^'''^]V6t;,.,«,0L2(T2) 
= 0(1) + oiRih)-^h^) = 0(1) + o(/ii+5-2"), 



which vanishes if we take a < . □ 



3+<5 
4 ■ 

Proof of Proposition 12.1. To prove Proposition 12.1, we first note that 



('^alT2x(A^\{0})xSA' WA)A1,(T-T2x§A)X«(r*T2xSA) ~ ^' 

since is ((?!)^)-invariant and (^^q , ^)^^(7..t2xSa)X''(t*t2xSa) ^ ^" ^hen, the (/)^-invariance of 
implies that i'q |T2x(Ai\{o})xSA vanishes. □ 
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Proof of Proposition 8.2. Proposition 12.1 implies that 

(At|T2xA^,a)^(j,.Tj2) .^o(T.T2) = 1™ (^OP/i [x {^—0^ a(a^.O^ Vh,v,^ 
for all a < and a £ "^^{^T^). The same holds if we replace x by : 



L2(T2) 



Since 



Op. (x^ (^) = Op, (x (^)) Op, (a) Op, (x (^)) (12.3) 

we obtain 

(//|T2xA-L,a)^(j,.T2),<g'0(T-T2) = 1™ (^Op,(a) Op, (^x (^^^)) ^/I'OPh (^X (^^^)) 



L2(T2) 



for aU Of < ^ and a £ '<f^{T*T^). □ 

13 Proof of Proposition 8.3: existence of the cutoff function 

Given a constant cq > 0, we define the following subsets of T^: 

= ({6 > co/i})a , ^, = / y B{x,{cohf')\ = [j B{x,{cohf'), '^h^-^h\^h, 

\xe{b>CQh} I ^ xeSh 

where for U C T^, we denote {U)^ := UreKi^ + '''"'} ^^"^ some cr e A-'- \ {0}. Remark that C 
and that = ^, U 5^, U (T^ \ Note also that the sets are non-empty for h small enough, 

and that is non empty (for h small enough) as soon as b vanishes somewhere on (this condition 
is assumed here since otherwise, GCC is satisfied). 

In this section, we construct the cutoff function Xh needed to prove the propagation results of 
Section 10. In particular, this function will be A-'--invariant and will satisfy Xh = on ^, and X;^ = 1 
on T2 \ 

The proof of Proposition 8.3 relies on three key lemmata. The first key lemma is a precised 
version of Proposition 6.1 concerning the localization in T*T^ of the semiclassical measure ^. It is 
an intermediate step towards the propagation result stated in Lemma 13.2. 

Lemma 13.1. For any x G ^c°°(M), .such that x = 1 « neighbourhood of the origin, for all 
a e '^^°^{T*t^), and 7 < ^^ , we have 



(Op,(a)w,,«;,)^2(T2) = (^Op;,(a)Op;, ^X (^^^^^-7^^^ ^ ^ +o(/i'^' '')|1 Op,(a)||£(L2), 

^ (13.1) 

For all a G '^^{T*T'^) and all r e R, 

{Oph{a o (j)r)wh,Wh)L2i^j2) = (Op,(a)w,,-u;,)^2(T2) +o(r/i^)|| Op,(a o </)i)|lioo(o,r;£(L2(T2))) 
In this statement, we used the notation 

II Op,(ao 0t)||ioo(o.^;£(L2(T2))) := sup \\ Opf^{a o (j)t)\\^(L^T^))- 

In turn, this lemma implies the following transport property. 
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Lemma 13.2. Suppose that the coefficients a,e satisfy 

< lOe < a, and a + 2e<l. (13.2) 

Then, for any time r G M uniformly bounded with respect to h, and any h-family of functions 
i/, -0^ e ^^°°(T2) satisfying 

lia^V'IU-(T^) < Cfc/i-2e|fe|^ forallk(^n\ (13.3) 

we have, 

{ip{s,y)wh,Wh)i^2(^j2-) = (?/'(s + T,y)wh,Wh)^2(T2) + (^'(s - t, y)w/i, Wh)i2(T2) 

+ O(/i"-i0^) + 0{h^-"-^') + o{h^), (13.4) 

where the coordinates {s,y) are the ones introduced in Section 7.3. 

In view of Proposition 8.3, this lemma will allow us to propagate the smallncss of the sequence 
u>/j above the set {b > Cq/i} to all S'f^. 

The third key lemma states a property of the damping function b, as a consequence of Assump- 
tion 2.12. 

Lemma 13.3. There exists bo = foo(£) > such that for all x G satisfying < b{x) < &o o,nd for 
all z g B{x, b{x)'^^), we have b{z) > 

With these three lemmata, we are now able to prove Proposition 8.3. 

Proof of Proposition 8.3. In the coordinates (s, y) of Section 7.3, we can write 

= TTa^ X Eh, = Ta^ X Fn, with E,, C F^ C Ta. 

Here, Eh is a union of intervals and has uniformly bounded total length. We can hence cover Eh with 
C\h~^' subsets of length of order (co/i)^'/2, overlapping on intervals of length of order (co/i)^'^/10. 
Associated to this covering, we denote by (V''j)j6{i,...,J}j = -^(K)-, ^ smooth partition of unity on 
Eh, satisfying moreover 

• E/=i^j(2/) = lforyGi;;,; 

• I|9^V',I|l-(t^) < Crnh-^'"\ for aU m £ N; 

• J = J{h) < Ch-^^. 

Similarly, we cover Ta-l with C2h~^^ subsets of length of order {cqH)'^' /2, overlapping on intervals 
of length of order (co/i)^'^/10, and define ('0fc)fce{i.. ..,/<■} a-n associated partition of unity on Ta^ 
satisfying 

• ^k{s) = 1 for s G Ta-l; 

• K ^ K{h) < Ch-^^, 

• for any k,ko & {1, . . . , K}^, there exists Tk satisfying \Tk\ < Lcngth(TA-L) < C and i^k{s + Tk) — 
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We set 

J K 

V'fe,(s,y) :=Vfc(s)^,(y), and xi^(s, y) = 1 - E E y) ^ "^°°(^')' 

]=i k=i 

which satisfies 9sX^(.s, y) ~ 0, i.e. Xh A^-invariant, together with 

• ~ on ff/i and hence h < coh on supp(x/^); 
. X;^ = 1 on \ 

• xt e [0, 1] on %, with \dyx^\ < Ch-^' and l^^x^^l < Ch^^'. 

To conclude the proof of Proposition 8.3, it remains to check Item 2 (|j(l — Xh)^h\\L^{T^) = o(l)), 
Item 4 {\\dyXhWh\\L^^T^) = o(l)) and Item 5 (Wd^Xh'^hWmT^) = o(l)). 

Now, let us fix jo G {1, ■ • ■ , Because of the definition of the set (^h, there exists fcp £ {1, . . . , A'} 
and xq G {b > cqH} such that supp('(/'fcojo) ^ B{xo, (coh)^^). According to Lemma 13.3, we have 
B{xo, {cohf") c {6 > ^}, so that supp(V'fcojo) C {6 > This yields 

^ii^kojoWh,Wh)L2{T^) < (&V'fcoJo'f«'i''^/i)i2(T2) = o(/l^+*), 

and hence {ipkojo^h, Wh)i^2(j2-^ ~ o{h^). Moreover, for any k e {1, . . . , K}, there exists T}^ satisfying 
|rfc| < 6*2 with 

Hence, using (13.4), we obtain 

0{h^) = {■llJko3ois,y)wh,Wh)L^T^) = {lpk3o{s + Tk,y)wh,Wh)L^T^) 
= ii^kjois + 2Tk,y)wh,Wh)L^T^} + {lpk3o{s,y)wh,Wh)L^T^) 

+ + 0(/ii""-2e) ^ o{h^). (13.5) 



Since both terms on the right hand-side are nonnegative, this implies {ipkjois, y)wh, w/i)l2(t2) = o{h^) 
as soon as 

a - lOe > 6, 
l-a-2e> S, 

(which implies (13.2)). From now on we will take S ~ 8e (this choice is explained in the following 
lines). The existence of a satisfying this condition together with (8.1) and a < 3/4, is equivalent to 
having e < 

To conclude the proof of Proposition 8.3, we first compute 

J K 

((1 - Xh)'>^h,Wh)L^(T2-^ = ^^{lpk]Wh,Wh)L^T2) = Ch~'^^o{h^) = oil), 

3 = 1 k=l 

since 5 > 4e. This proves Item 2. Next, we have by construction supp(9^Xh) ^^w{dyXh) C 
'^h with ||5,x;^I|l~(t=) = 0{h-^^), Wdlx^L^ij^) = 0{h-^n. Hence, covering supp(9,xi^)) by 
balls of radius (cq/i)"" and using a propagation argument similar to (13.5) shows that we have 

l|w''i|lL2(supp(a„x;^)) = °(^^)- We thus obtain 

lia,xN/.IU^(T^) = oih^-'^) = o(l), \\dlx^iw,h2^j2^ ^ o{h^-^n = 0(1), 
(since 5 > 8e) which concludes the proof of Items 4 and 5, and that of Proposition 8.3. □ 
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To conclude this section, it remains to prove Lemmata 13.2, 13.1 and 13.3. In the foUowing 
proofs, we shall systematically write 77 in place of P\(^ and cr in place of (1 — Pa)^ to lighten the 
notation. Hence, ^ g is decomposed as ^ = rj + a with G (A) and cr e A-'-, in accordance to 
Section 7.3. 

Proof of Lemma 13.2 from Lemma 13.1. First, given a function ip £ ^^°°(T^) satisfying (13.3), we 
have, 

(V'W/i, W/i)i2(T2) = (Op^(V' O (f>r)wh, W/i)i2(T2) + 0{Th^)\\ Opi^ii' O 0t)IU°°(O.r;£(L2)) 



Opft(V'°0r)Op,, X 



hi 



L2(T2) 



+ {0{Th^2')j^0{Th''^2' 7)) II Op;,(V'0 0t)IU~(O,r;£(L2)), 



when using Lemma 13.1 together with Op^ (x (2^^)) '^h ~ Wh- Next, the pscudodifFcrential calculus 
yields 



3 + 3 



_n_ 



L2(T2) 



+ (o(t/i 2 )+o(r/i 2 '^))|| Op^(V^o ,/,t)||ioo(o^^;£(L2)). 



(13.6) 



A particular feature of the Weyl quantization in the Euclidean setting is that the Egorov theorem 
provides an exact formula (see for instance [DS99]): Op,j('0 o cp^) = e^'*'*" Opf^{ip)e^*^~ , so that 
II ^Phi''P ° 4't)\\ L°° {o,t;C(l^)} — ^0 uniformly with respect to h. Now, remark that the cutoff function 



-) X ('TT^) '^^^ decomposed (for h small enough) as 



X 



V 



X 



for some nonnegative function xjj such that (cr, rj) 1-^ xl^if^) G '^^^{^^), such that xj;(''') = X ('^^W^) 

for 77 G suppx (2^) ^nd cr > 0, and x5j(cr) = for 77 ^ suppx (271^) oi" cr < 0. 
Choosing 7 = a, we have in particular 

\a-l\<Ch" on supp(x(^)x:;(^))- 

Next, we recall that ipo(j).,.[s, y, cr, 77) = iJj{s + Ta, y + Trf), and we focus on the first term (corresponding 
to cr > 0) in the right-hand side of the identity 



(13.7) 



We set 



C|^^(s,2/,cr,7/) = X (^) x!5(o-)V'(s + rcr,y + r77), and Cr^^ (s, y, cr, 77) = x (^) xJ,'(o-)V'(s + "t, y), 
and wc want to compare Op/j(Ct^'') and Op^(Cl^''). For this, let us estimate, for multiindices £,m £ N^, 

<CmY. 1^:.) (^) X:;(^)) ^(...)^(<.,.) (^(^ + 2/ + rv) i^{s + r, y))| . (13.8) 



On the one hand, we have 



(13.9) 
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On the other hand, for \v\ > we can also write 



9(s,y)9(a,rj) + T<T,y+ TTj) - V(s + T, J/)) = dl^ y)d^„ ,^-^^p{s + TCT,y + T7j) 

< C,,,|r|l''l/l-2-(l^l + l''l) < + 

since |r| < C. 

Finally, for |:/| = 0, we apply the mean value theorem to the function 



(cr, 77) h-> df^ y^ij{s + Ta,y + T-q) 



and write 



^(s,y) i^i^ + T(J,y + Trf) ~ 'ijj{s + T, y)) 

< (hi + |ct - 1|) sup V(^,^)af )(V'(s + rcr,y + r7/)) 

With (13.3), this yields 

dts,y) + ra,y + ttj) - ^{s + r, y)) | < M + k - l|)Q/j-2e|^| |^|;^- 

< (|r;| + |a-l|)Q/i-2-(l^l+i), 



(13.10) 



for |t| < C. 

Using now that I??] < C/i" and |cr-l| < C/i" onsupp (x (5^^) x^J(ct)) , and combining (13.8), (13.9) 
and (13.10), we obtain, for aU m £ e and < h < ho sufficiently small, 



0<i^<m 



2e 



Using a precised version of the Calderon-Vaillancourt theorem, as presented in Theorem A.l below 
(in which only \£\ = 4 derivations are needed with respect to x in dimension two), we obtain 



Op,(c(^)) = Op,(e)) + o^(i.)(/."-i"^). 



Similarly, we have 



Op, [x[t^) x';{-<Tms + ra,y + Tq)) = Op, [x [-^] x'^h'^ms - r,y)] + OciL^^ih^-'^^. 



2h" 



2h" 



Coming back to (13.6) and using (13.7), we finally obtain, for all |t| < C, 

i4>Wh, w/i)l2(t2) = (Op, (^) xJj(c^)V'(s + T, y)j Wh,Wh^^^ 



+ (Op. (x (^) X^(--)^(^ - r,y)) -/u-.)^^ 



(T^) 

+ ©(/l"^!"^) + ©(/l^-^-^e) _^ o{h^) + 0{h^-°'). 

With the pseudodifferential calculus, this yields (13.4), which concludes the proof of Lemma 13.2. □ 

Proof of Lemma 13.1. Here, we only have to make more precise some arguments in the proof of 
Lemma 6.1. Recall that according to Lemma 8.1, Wh satisfies Pj^Wh ~ o{h?^^). 
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First, we take x G '^^{^), such that x = 1 in a neighbourhood of the origin. Hence, "'" ^^^^ € 
"(R) and we have the exact composition formula 



since both operators are Fourier muhiphers. Moreover, Op^ — x \ — j j uniformly 
bounded as an operator of £(L^(T^)). As a consequence, we have 



Op/i(a) Op,, ( 1 - X ( ^^^^-7^ ) ) Wh,Wh 



L2(T2) 



Op.(a) Op, ((1 - X (^)) ^) 



L2(T2) 



where A = Op, (a) Op, ((l - X (^^^)) IcF^) bounded on ^2(^2). Using Pl'wh = o{h:^+^) 
and w,)i2(T2) = o{h^+^), this gives 



(^Op,(a) Op, (^1 - X (^^)) '"'''^'O^,^^,^ ^ Op,(a)|U(i2), 
which in turn implies (13.1). 



Next, Identity (6.2) yields, for aU a G '^^^(T'^), 

(Op,(C ■ d^a)wh,'Wh) L2(^j2) = ^ {Opf^{a)wh, PbWh) ^2i^j2-i - ^ (Op,(a)Pb''w,, u;,)^^^.^,^^ 

^ ^ (Op,(a)u;,,6u;,)^2(T2) - ^ (Op,(a)5u;,, w;,)^2(t2) 



o(/i^+'>)|| Op,(a)||£(i2) + o(/i — )|| Op,(a)||£(i2), 

as a consequence of P/^Wh = o{h'^~^^) and (few,, z/;,)^2(t2) = o(/i^^*). Applying this identity to a o i/ij 
in place of a, and integrating on i € [0, r] finally gives 

(Op,(ao (^^)w,,u;,)^2(t2) = (Op,(a)w,,w,)^2(T2) +o(r/i^)|| Op,(a o (/it)||ioo(o,r;£(L2)), 

which concludes the proof of Lemma 13.1. □ 

Proof of Lemma 13.3. Here, B B{x,b{x)'^^) denotes the euclidian ball in centered at x of 
radius b{x)^^ . Setting 

M := sup6(z), m := inf b(z), 

we have 

|V6(z)| < C^b^-^z) < CeM^-^, for all z e B, 
as a consequence of Assumption 2.12. Moreover, the mean value theorem yields 

bix) - CeM^-^bixf^ < b{z) < b{x) + C'eM^-^bixf^ for z e B, 
and, in particular, 

m>b{x)-CeM^-'b{xf', and M < b{x) + CeM^^'bixf" . (13.11) 
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Now, defining f{M) := b{x) + CgM^ we see that / is a strictly concave function with 

/(O) = b{x) > 0. There exists a unique Mq £ R+ satisfying /(A/q) = M^. Moreover, we have 
M < f{M) if and only if M < Mq. Taking bo sufficiently smah so that 60 + C'e^o"*""' ^ ^a'' ^ 
we obtain /(6(a;)^~'^) < b{xY^^ . In particular, this gives Mq < b{xy~^ and hence M < b{xY^^ 
according to the second estimate of (13.11). Coming back to the first estimate of (13.11), this yields 

TO > b{x) - C,b{xY^~'^'Hx?' = b{^) ~ C.bi^y^"'- 

Taking bg sufficiently small so that bo — C^bg'^^ > we obtain m > which concludes the proof 
of Lemma 13.3. □ 



Part IV 

An a priori lower bound for decay rates 
on the torus: proof of Theorem 2.5 

Under the assumption 

{6>O}n{a;o + reo,re]R}==0, (13.12) 

for some {xo,£,o) € T*T^, 7^ 0, we construct in this section a constant kq > and a sequence 
{(Pn)neN of C'(l)-quasimodes in the limit n — >■ +oo for the family of operators P^inno). 

We use the notation introduced in Sections 7.1 and 9. First, note that, as a consequence of (13.12), 
^0 is necessarily a rational direction, and the set {a;o + t^q,t € K} is a one-dimensional subtorus of 
T^, given by 

{a;o + T^o, T e M} = {xo + t^o,t e M} = xq + Tj^± , with e V. 

Let X G '^^c^(T^) such that x ha-s only x- Fourier modes in A^„, x = on a neighbourhood of 
{b > 0} and x = 1 on .tq + Tj^± . 

So 

From Assumption (13.12), we have rk(A^g) = 1, so that one can find k e Aj-^ r\I?\ {0}. Besides, 
for all n G N we have nk e A^^ n \ {0}. 

Wc then define the sequence of quasimodes {^n)n&i by 

^„(x) = x(a;)e"'=•^ n e N, x eT\ 

We have e "^""(T^). together with the decoupling 

(p„ o ^A^^ (s, y) = x(2/)e™'=•^ n e N, (s, y) £ T^^^^ x Ta,„ . 

This yields 

- (JAj„ ATXj^ ) O TTAj^ (S, ?/) = - (Aa^^ + AaJ-^ ) O TTA^^ (s, y) 

= -e™'=-AA,„x(y) + r^2|fc|\(2/)e™'=•^ 

Moerover, 6(p„ = 0, according to their respective supports. Hence, recalling that P{in\k\) = —A — 
n^jfcp + in\k\b{x), we have 

(Ta,^ P(^n|fc|)TX^J^„ o tta,^ = -e™'=-AA,„ x(y), 

and 

||P(m|fc|)(p„||i2(T2) = \\{Tf,^^^P{in\k\)T;{^Jip„ o TTA^JmT^^ xTa^J = CoH Aa^^xIU^ITa^J- 
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Since we also have ||<^„||i2(T2) = |1Ta (y9„j|i2(T , xTa, ) = C'o|1xI1l2(Ta, ), we obtain, for all n e N, 



\\Wn\\j2n2\ X L2(Ta, ) 

which concludes the proof of Theorem 2.5. □ 
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A PseudodifFerential calculus 

In the main part of the article, we use the semiclassical Weyl quantization, that associates to a 
function a on r*R^ an operator Op/j(a) defined by 

(Op.(a).)(.) I e^^ (^-)a uiy)dy d^. (A.l) 

For smooth functions a with uniformly bounded derivatives, Op/j(a) defines a continuous operator 
on ^(M^), and also by duality on ^'(M^). On a manifold, the quantization Op^ may be defined by 
working in local coordinates with a partition of unity. On the torus, formula (A.l) still makes sense : 
taking a € '^°°(T*T^) is equivalent to taking a € "^"^ (M? x R^) , (27rZ)2-periodic with respect to the 
x-variable. Then the operator defined by (A.l) preserves the space of (27rZ)^-periodic distributions 
on and hence ^'(T^). 

We sometimes write, with D := id, 

a{x, hD) ~ Op^(a). 

We also note that Op^ (a) is the classical Weyl quantization, and that we have the relation 

a{x,hB) = OpJa(.T,0) = OvxH^M))- 

Theorem A.l. There exists a constant C > such that for any a G '^°°(T*T^) with uniformly 
bounded derivatives, we have 

II Opi(a)||£(i2(T2)) < C ||a^a^a||ioo(r.T2). 

a,/3G{0,l,2}2 

Equivalently, this can be rewritten as 

l|0p,(a)|U(i2(T.)) <C h\^\\\dy^^aU^^r'T^y 
Q,,3e{04,2}2 

This precised version of the Calderon-Vaillancourt theorem is needed in Section 13, and proved 
in [Cor75, Theorem Bp] or [CM78, Theoreme 3]. Here in dimension two, this means that only 
|a| = 4 derivations are needed with respect to the space variable x. 

B Spectrum of P{z) for a piecewise constant damping 
(by Stephane Nonnenmacher) 

In this Appendix we provide an explicit description of some part of the spectrum of the damped wave 
equation (1.1) on T^, for a damping function proportional to the characteristic function of a vertical 
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strip. We identify the torus with the square {—1/2 < x < 1/2, < y < 1}. We choose some 
half- width tr € (0, 1/2), and consider a vertical strip of width 2a. Due to translation symmetry of 
T^, wc may center this strip on the axis {x — 0}. Choosing a damping strength B > 0, we then get 
the damping function 

bix,y)=bix) = ll 't^"' (B.l) 
\B, cr < |a;| < 1/2 . 

The reason for centering the strip at .t = is the parity of the problem w.r.t. that axis, which greatly 
simplifies the computations. 

We are interested in the spectrum of the operator A generating the equation (1.1), which amounts 
to solving the eigenvalue problem 

P{z)u = 0, for P{z) = -A + zb{x) + z^, z e C , u e L^{T^) , u ^ . (B.2) 

This spectrum consists in a discrete set {zj}, which is symmetric w.r.t. the horizontal axis: indeed, 
any solution (z, u) admits a "sister" solution (z, u). Furthermore, any solution with Imz 7^ satisfies 

1 (u, 6w)r2CT2'l ~ 

Rez = — ^ ^ ^ and thus -B/2<Rcz<0. (B.3) 

We may thus restrict ourselves to the half-strip {—B/2 < Rcz < 0, Imz > 0}. 

Our aim is to find high frequency eigenvalues (Imz 3> 1) which arc as close as possible to the 
imaginary axis. We will prove the following 

Proposition B.l. There exists Co > such that the spectrum (B.3) for the damping function (B.l) 
contains an infinite subsequence {zi} such that Imz^ — )■ 00 and \ He Zi\ < ji~~^3j2 ■ 

The proof of the proposition will actually give an explicit value for Co, as a function of B, a. 
Proof. To study the high frequency limit Im z — >■ 00 we will change of variables and take 

z^i{\lh-vl), 

where h <E (0, 1] will be a small parameter, while C <= C is assumed to be uniformly bounded when 
h ^ The eigenvalue equation then takes the form 

(-/i^A + ih{\ + hC)h)u = (1 + 2/iC(l + hZl2))u. (B.4) 

Having chosen h independent of y, wc may naturally Fourier transform along this direction, that is 
look for solutions of the form u{x,y) = e^'^'^^'^ v{x) , n G Z. For each n, we now have to solve the 
1-dimensional problem 

{~h^d^/dl + ih{l + hC,)h{x))v = (1 - {2TThnf + 2K(1 + hC/2)) v . (B.5) 

Let us call 

B='B(l + /iC), C C(i + /iC/2) . 
In terms of these parameters, the above equation reads: 

{~h^d^/dl+ihBl[^^l^l<i/2}{x))v ^ Ev , with E = 1 - {2TThn)^ + 2hC . (B.6) 
Since we will assume throughout that ^ = 0(1), we will have in the semiclassical limit 

B = B + Oih), C = C(l-K/2 + 0(/i')). (B.7) 

At leading order we may forget that the variables B, C are not independent from one another, and 
consider (B.6) as a bona fide linear eigenvalue problem. 
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Since the function b{x) is even, we may separately search for even, resp. odd solutions v{x). Let 
us start with the even solutions. Since b(x) is piecewise constant, any even and periodic solution 
v{x) takes the following form on [—1/2, 1/2] (up to a global normalization factor): 

^ fcos(fca;), |x| < cr, 

"^""^ \/3cos(fc'(l/2- a<\x\<l/2,' 

E^/^ (E-ihBV/^ 
k = ^, k' = ^ ' . (B.9) 

h h 

Wc notice that k,k' are defined modulo a change of sign, so we may always assume that Refc > 0, 
Refc' > 0. The factor (3 is obtained by imposing the continuity of v and of its derivative v' at the 
discontinuity point x ~ a (we use the notation cr' 1/2 — cr): 

cos{ka) = /3cos(fc'cr'), 
-k sm{ka) = 13k' sin(fc'cr') . 

The ratio of these two equations provides the quantization condition for the even solutions: 

k' 

tan(fccr) = tan(fc'cr') . (B.IO) 

k 



Similarly, any odd eigenfunction takes the form (modulo a global normalization factor): 
v{x) 



sin(fc a;), |a;| < cr, 

(3agn{x) sin(fc'(l/2- |x|)), cr < \x\ < 1/2, 



so the associated eigenvalues should satisfy the condition 

k 

tan(/ecr) = -— tan(fc'CT') . (B.ll) 
k 

We will now study the solutions of the quantization conditions (B.IO) and (B.ll), taking into account 
the relations (B.9) between the wavevectors k, k' and the energy E. To describe the full spectrum 
(which we plan to present in a separate publication), we would need to consider several regimes, 
depending on the relative scales of E and h. However, since we are only interested here in proving 
Proposition B.l, we will focus on the regime leading to the smallest possible values of | lm(\ = | Rez|. 
What characterizes the corresponding eigenmodes v{x) 1 From (B.3) we see that the mass of v{x) in 

1 /2 

the damped region, 2 \v{x)\'^ dx, should be small compared to its full mass. Intuitively, if such a 
mode were carrying a large horizontal "momentum" Re(/ifc) in the undamped region, it would then 
strongly penetrate the damped region, because the boundary at a; = cr is not reflecting. As a result, 
the mass in the damped region would be of the same order of magnitude as the one in the undamped 
one. This hand-waving argument explains why we choose to investigate the eigenmodes for which hk 
is the smallest possible, namely of order 0{h). This implies that E = [hkY = 0{h^), which means 
that almost all of the energy is carried by the vertical momentum: 

hn = {2Tr)-^ +0{h) . 

The study of the full spectrum actually confirms that the smallest values of Im C are obtained in this 
regime. 

Eq.(B.9) implies that the wavevector k' in the damped region is then much larger than k: 
k' . i-^'^B+jhkn/^ ^ e-/\B/hy^^ + 0(/.V^) . 

Imfc'cr' « -cr'(B/2/i)i/2 is negative and large, so that tan(fc'cr') = + ©(e^^^^'''^')), uniformly 
w.r.t. Re(fc'cr'). 
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Even eigenmodes 

In this situation the even quantization condition (B.IO) reads 

tan(fca) = ^^ {l + ©(g-^'^^^/'')''')) . (B.12) 

Since the r.h.s. is large, ka must be close to a pole of the tangent function. Hence, for each integer 
m in a bounded interval^ < to < M we look for a solution of the form 

7r(TO + l/2) „ , ^ 

«m+l/2 = h dfcm+1/2, With |()fc„+i/2| < 1 . 

(T 

The quantization condition (B.12) then reads 

Using (B.6), the corresponding spectral parameter C is then given by 



Cn,m+l/2 



2h 

{'^^hnf - 1 /^(to + 1/2) X 2 /2 /7r(TO + l/2)x 2 

2/1 + 2I a J + V a ) +^1^^ ^ 



From the assumptions on the quantum numbers ri, to, we check that Qn.m+1/2 ~ ^{^)- We may now 

go back to the original variables C, B, using the relations (B.7). The spectral parameter C has an 
imaginary part 

ImCn,™+l/2 = ImC„,™+l/2(l - /lReC„,™+l/2) + Oih") = /,3/2 (^("^ + V^))' ^ 0(/^2) ^ (g^^3) 

cr'^(2i3)^/^ 

Returning back to the spectral variable z, the above expression gives a string of eigenvalues {zn^m+1/2} 
with Imz„ „i+i/2 = /i~^ + 0{1), Rez„ „j+i/2 = ~ Im Cn,m+i/2- These even-parity eigenvalues prove 
Proposition B.l, and one can take for Co any value greater than \ . □ 

We remark that the leading order of fcm+1/2 corresponds to the even spectrum of the operator 
—h?d'^/d'^ on the undamped interval [— ct, tr], with Dirichlet boundary conditions. The eigenmode 
Vn m+i/2 associated with Cn,m+i/2 is indeed essentially supported on that interval, where it resembles 
the Dirichlet eigenmode cos (a;7r(l/2 + to)/(t). At the boundary of that interval, it takes the value 

Vn,m+U2{<^) = (-l)™+ie^3-/4 j^y, 7:{m + l/2) ^ ^^^^ 

and decays exponentially fast inside the damping region, with a "penetration length" (Im/c')^^ « 
{2h/Bf/'^. From (B.3) we see that the intensity |'yn,m+i/2(o')P ^ Ch penetrating on a distance 
^ h^/"^ exactly accounts for the size ^ h^^"^ = hh^^^ of the Re Zn^jn+i/2- 

We notice that the smallest damping occurs for the state Vn,i/2 resembling the ground state of 
the Dirichlet Laplacian. 



^Recall that we only need to study values Rcfc > 0. 
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Odd eigenmodes 

For completeness we also investigate the odd-parity eigenmodes with k = 0(1). The computations 
are very similar as in the even-parity case. The odd quantization condition reads in this regime 

tan(fca) = »^ (l + 0{e-^'-^ . (B.14) 

The r.h.s. is then very small, showing that ak is close to a zero of the tangent, so we may take 
km = nm/a + 5km with \5km\ ^ 1 and < m < M . We easily see that the case m = docs not 
lead to a solution. For the case m > we get 

and thus 

^ra = —\\^-h^l^—^^+0{h)y \<m<M. 

These values fc^, approximately sit on the same "line" {s(l + h^l"^ ^^rp^)-, s G M} as the values 
km+i/2 corresponding to the even eigenmodes, both types of eigenvalues appearing successively. The 
corresponding energy parameter Cn,m satisfies 

ImC„.^ = /.3/^ ir5u/2 +^(^')- (^-15) 

As in the even parity case, the eigenmodes Vn,m are close to the odd eigenmodes sin (^xmn/a) of the 
semiclassical Dirichlct Laplacian on [—a, a], and penetrate on a length ^ ft,^/^ inside the damped 
region. 

The case of the square 

If the torus is replaced by the square [—1/2,1/2] x [0,1] with Dirichlet boundary conditions, with 
the same damping function (B.l), the eigenmodes P{z) can as well be factorized into u{x,y) = 
sv[\{2TTny)v{x), with n S \ 0, and v{x) must be an eigenmode of the operator (B.6) vanishing at 
X = ±1/2. We notice that the odd-parity eigenstates (B) satisfy this boundary conditions, so the 
eigenvalues z^.m (with real parts given by (B.15)) belong to the spectrum of the damped Dirichlet 
problem. 

Similarly, in the case of Neumann boundary conditions the eigenmodes factorize as u(x, y) = 
cos{2Tmy)v{x), with n G ^N. The even-parity states (B.8) satisfy the Neumann boundary conditions 
at a; = ±1/2, so that the eigenvalues described in (B.13) belong to the Neumann spectrum. 

As a result, the Dirichlet and Neumann spectra also satisfy Proposition B.l. 
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